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Abstract  of  Dissertation  Presented  to  the 
Graduate  Council  of  the  University  of  Florida  in  Partial 
Fulfillment  of  the  Requirements  for  the  Degree  of  Doctor  of  Philosophy 

STABILITY  OF  HETEROGENEOUS  SWIRLING  FLOWS 

By 

Yee-tak  Fung 
December,  1974 

Chairman:  Dr.  Ulrich  H.  Kurzweg 

Major  Department:  Engineering  Sciences 

A general  mathematical  investigation  of  the  hydro- 
dynamic  stability  of  heterogeneous  swirling  flows  with 
possible  applications  to  the  gaseous  centrifuge  method 
for  separation  of  uranium  isotopes  and  to  the  proposed 
vortex-driven  gas-core  nuclear  reactor  is  presented.  The 
flow  under  consideration  is  assumed  to  be  inviscid  and  in- 
compressible but  has  radially  dependent  density  and  velocity 
components  within  an  annular  region  between  two  concentric 
cylinders.  Gravitational  effects  are  neglected.  The  con- 
cept of  "centrifugal  stability"  is  introduced  to  distin- 
guish instabilities  caused  by  centrifugal  forces  due  to 
rotation  of  the  flow  from  those  caused  by  shearing  effects 
due  to  the  relative  motion  between  two  adjacent  fluid  layers. 
A flow  distribution  is  "centri fugal ly  stable"  if  the  product 
of  its  density  and  the  square  of  its  circulation  is  an 
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increasing  function  of  radius.  It  should  be  noted  that 
the  classical  Rayleigh- Synge  criterion  applies  only  to 
centri fuga 1 1 y induced  instabilities.  A generalized  suf- 
ficiency condition  for  hydrodynamic  stability  of  swirling 
flows  against  arbitrary  infinitesimal  perturbations  is  ob- 
tained. The  sufficiency  condition  states  that  a non-di- 
mensional parameter  A,  reminiscent  of  the  familiar  Richard- 
son number  encountered  in  two-dimensi ona 1 stratified  shear 
flows,  shall  be  greater  or  equal  to  a quarter  at  all  points 
within  the  flow  for  stability. 

Extensions  of  Howard's  semi-circle  theorem  are  ob- 
tained for  four  types  of  swirling  flows  with  centri fugally 
stable  velocity  and  density  distributions.  The  results 
are  compared  with  the  exact  solutions  obtained  for  three 
special  velocity  and  density  profiles.  A semi-circle  bound 
for  instability,  growth  rates  and  phase  velocities  is  also 
obtained  for  solid  body  rotation  with  arbitrary  density 
distribution  when  the  flow  is  subject  to  azimuthally  peri- 
odic disturbances. 

Two  interfacial  conditions  for  di  sconti  nuous.  densi  ty 
and  velocity  distributions  are  obtained  via  both  mathematical 
and  physical  arguments.  In  particular,  the  centrifugal  ef- 
fects due  to  discontinuities  cannot  be  neglected  when 
matching  techniques  are  applied  at  interfaces.  Effects  of 
surface  tension  on  stability  characteristics,  which  are 
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found  to  always  stabilize  the  flows  except  for  axisym- 
metric  modes  with  axial  wave  numbers  less  than  a certain 
value,  become  important  for  flows  with  one  or  more  inter- 
faces. 

Four  types  of  velocity  and  density  profiles  are 
examined  analytically,  and  exact  solutions  expressed  in 
terms  of  the  modified  Bessel  functions  are  obtained.  Sta- 
bility characteristics  of  these  flows  are  found  and  shown 
to  always  satisfy  the  general  criteria  derived  earlier. 

It  is  found  that  pure  rotational  flows  without  axial  veloc- 
i ty  gradients  with  large  shear  are  more  unstable  to  two- 
dimensional  azimuthally  periodic  disturbances  than  either 
axisymmetric  or  more  general  three-dimensional  disturbances. 
Furthermore,  solutions  with  large  wave  numbers  are  most  un- 
stable for  velocity  profiles  containing  one  or  more  vortex 
sheets. 


Chai rman 


XV 


CHAPTER  1 


INTRODUCTION 

I . 1 Physical  Considerations 

A recent  report  by  Abajian  and  Fishman  (1973)  on 

the  supply  of  enriched  uranium  indicates  the  possible 

importance  of  the  high-speed  gaseous  centrifuge  process 

235 

for  separation  of  uranium  isotopes.  Isotope  U is  sep- 

2 3 8 

arated  from  isotope  U in  a high-speed  rotating  cylinder 
because  of  the  slight  difference  in  molecular  weights. 

The  gas  molecules  experience  outward  centrifugal  forces 
produced  by  the  rotation,  and  are  subject  to  thermal  dif- 
fusive forces  which  tend  to  keep  the  molecules  evenly  dis- 
tributed. Since  the  diffusion  effect  is  smaller  for  the 
heavier  molecules  there  are  rel ati vely  fewer  specimens  of 
the  heavier  isotope  at  the  inner  radius.  A major,  and  as 
yet  incompletely  considered  problem,  associated  with  the 
centrifugal  separation  process  concerns  the  hydrodynamic 
stability  of  this  flow.  For  isotope  separation  it  is 
clearly  a requirement  that  the  flow  be  stable  for  all  pos- 
sible hydrodynamic  disturbances. 

The  need  for  a stability  analysis  on  a rotating 
flow  of  variable  density  is  also  seen  in  the  vortex-driven 
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gas-core  nuclear  reactor  proposed  by  the  Oak  Ridge  National 
Laboratory.  In  this  method  for  generating  nuclear  power, 
a gas  of  low  molecular  weight  is  injected  tangentially  into 
a cylindrical  chamber  along  with  fissionable  materials  also 
in  gaseous  form.  The  gas  is  withdrawn  either  through  a 
central  porous  tube,  or  through  the  open  ends  of  a central 
cylinder,  thus  generating  a potential  vortex  inside  the 
chamber.  The  gaseous  fuel  being  of  heavier  molecular  weight 
moves  radially  outward,  while  the  light  gas  swirls  radially 
inward  and  moves  toward  the  exit.  An  annular  cloud  of  gas- 
eous fuel  is  then  accumulated  at  some  radial  position,  and 
will  become  heat  producing  when  the  neutron  flux  of  the 
cloud  achieves  a certain  critical  value.  Again  this  is  a 
flow  situation  which  requires  careful  consideration  for 
the  condition  of  stability. 

The  recent  speculations  by  Soviet  and  American  re- 
searchers on  the  possibility  of  achieving  a fusion  reaction 
with  high-powered  lasers  has  led  to  several  proposals  for 
controlling  such  reactions.  One  such  design  for  a laser 
fusion  reactor  involves  containment  of  a laser  imploded 
hydrogen  pellet  with  a liquid  lithium  vortex,  the  liquid 
lithium  being  fed  tangentially  into  a spherical  pressure 
vessel.  A free-standing  vortex  is  formed  in  the  center 
as  a result  of  the  high  velocity  spiralling  lithium.  Dur- 
ing the  fusion  process,  the  energy  released  by  the  fused 
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deuterium  is  absorbed  in  this  liquid  lithium  vortex.  Due 
to  the  properties  of  metal,  the  lithium  near  the  center 
of  the  vortex  where  the  fusion  takes  place  has  a higher 
temperature  and  a lower  density.  The  hot  lithium  is  thus 
kept  in  the  center  of  the  spherical  pressure  vessel  by  es- 
sentially centrifugal  effects.  As  a result  the  heat  pro- 
duced can  be  transferred  by  extracting  the  high-temperature 
liquid  lithium  through  the  bottom  of  the  vessel.  This  ex- 
ample of  a heterogeneous  swirling  flow  is  different  physi- 
cally from  the  previous  two,  but  they  all  show  very  similar 
hydrodynamic  conditions  which  need  to  be  examined  for  their 
hydrodynamic  stability  characteristics.  We  would  like  to 
carry  out  such  a stability  analysis  for  heterogeneous  swirl- 
ing flows. 

1 . 2 Mathematical  Model 

The  study  of  the  hydrodynamic  stability  for  hetero- 
geneous incompressible  swirling  flows  with  radially  depen- 
dent angular  and  axial  velocity  components  and  density 
Pq('")  is  of  considerable  practical  interest  in  connection 
with  the  just-mentioned  physical  problems.  The  swirling 
stability  of  such  flows  is  characterized  by  the  exchange 
of  position  between  two  fluid  elements  in  a centrifugal 
force  field,  i.e.,  the  transfer  of  energy  from  the  mean 
velocity  field  to  small  amplitude  perturbations  in  the 
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form  of  travelling  waves.  The  instability  mechanism  is 
similar  to  that  characterizing  the  stability  of  hetero- 
geneous two-dimensional  parallel  shear  flows.  In  such 
shear  flows  the  stability  is  governed  by  the  magnitude  of 
the  ratio  of  the  density  gradient  to  the  shearing  rate, 
the  Richardson  number.  To  understand  the  behavior  of  this 
instability  mechanism  in  swirling  flows,  it  is  necessary 
to  consider  the  mathematics  of  the  Navi er-Stokes  equations, 
the  continuity  equation,  and  the  condition  of  incompressi- 
bility and  solve  these  governing  differential  equations. 

The  onset  of  instability  can  be  determined  by  assum- 
ing temporal  or  spatial  disturbances  whose  behavior  is 
governed  by  the  linearized  forms  of  the  governing  non-linear 
equations.  However,  one  should  keep  in  mind  that  instabil- 
ity does  not  necessarily  transfer  a laminar  flow  into  turbu- 
lence. Unstable  modes  sometimes  are  simply  a transition 
from  one  type  of  laminar  motion  to  another.  Discussion  of 
the  stability  beyond  the  initial  onset  is  not  within  the 
capability  of  a linear  analysis  and  will  not  be  considered 
here . 

It  has  been  accepted  that  the  method  of  normal 
modes  is  complete  and  well  explains  the  linear  wave  behav- 
ior in  fluid  mechanics  (see  Lin  1955,  Roberts  1967,  Betchov 
and  Criminale  1967  ).  The  linearized  partial  differential 
equations  under  consideration  are  cyclic  with  respect  to 
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time  and  to  two  spatial  variables.  This  strongly  suggests 
they  be  solved  by  Laplace  (time)  and  Fourier  (space)  trans- 
form methods  (Case  1960  a,  b).  The  differential  equations 
after  transformation  are  ordinary  and  in  general  complex 
both  with  respect  to  the  temporal  and  spatial  parameters 
(Gaster  1963).  If  any  one  of  the  Laplace  or  Fourier  modes 
is  unstable,  we  consider  the  system  to  be  unstable  as  a 
whole.  In  other  words,  this  method  is  assuming  that  each 
perturbation  can  be  solved  into  dynamically  independent 
wave  components.  In  our  analysis  we  will  restrict  our- 
selves only  to  temporal  perturbations. 

In  the  case  of  an  inviscid  analysis  as  will  be 
carried  out  here,  the  eigenvalue  problem  will  become  sing- 
ular, and  the  singularity  admits  solutions  with  discontinu- 
ous derivatives  and  a continuous  spectrum  of  eigenvalues 
in  addition  to  well-behaved  solutions  with  a discrete  spec- 
trum (Lamb  1932,  §151).  This  singularity  may  be  regarded 
as  a consequence  of  either  our  assumption  of  periodic  wave 
motions  or  our  neglect  of  viscous  (or  diffusion)  effects. 
Two  difficulties  arise  as  a result  of  such  an  analysis: 
first,  a critical  layer  is  introduced  inside  the  flow  do- 
main, which  can  be  smoothed  out  only  by  addition  of  viscos- 
ity (or  diffusion);  second,  the  solution  to  the  governing 
equations  has  as  a complex  eigen  wave  speed,  c = c^  + i c^ , 
and  its  complex  conjugate.  That  is,  both  damped  and  growing 
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waves  exist  simultaneously  for  the  same  wave  speed  and  it 
accordingly  becomes  much  more  difficult  to  locate  the 
stability  boundary  (refer  to  Chapter  2 for  definition) 
except  for  some  special  profiles  in  which  the  principle 
of  exchange  stabilities  is  valid  (Howard  1963). 

1 . 3 Previous  Work 

It  was  early  shown  by  Rayleigh  (1916)  that  the  in- 
viscid  stability  of  fluid  motion  in  cylindrical  strata 
under  axisymmetric  disturbances  required  only  that  the 
square  of  the  circulation  increase  outwards.  His  argument 
was  based  on  an  analogy  between  the  motion  and  the  equil- 
ibrium of  a heterogeneous  liquid  under  gravity  and  was  ac- 
complished by  the  consideration  of  the  increment  in  kinetic 
energy  consequent  on  an  interchange  of  two  rings  of  liquid. 
Under  this  criterion,  a homogeneous  potential  vortex  with- 
out axial  flow  gradient  is  neutrally  stable  against  axisym- 
metric disturbances.  Von  Karman  (1934)  gave  a slightly 
different  physical  derivation  of  the  Rayleigh  criterion 
by  letting  the  centrifugal  force  balance  the  pressure 
gradient  established  by  the  velocity  field. 

Synge  (1933)  gave  a rigorous  proof  of  the  Rayleigh 
criterion  by  taking  into  account  the  radially  dependent 
density  variation  and  casting  the  problem  into  the  form  of  a 
Sturm-Li ouv i 1 1 e system.  The  well-known  Ray 1 e i g h -Synge 


7 


criterion  then  states  that  an  inviscid  rotating  flow  is 
stable  against  axi symmetri c disturbances  whenever  the 
product  of  its  density  and  the  square  of  its  circulation 
is  an  increasing  function  of  radius. 

Ponstein  (1959)  first  looked  at  the  stability  of 
rotating  cylindrical  jets  for  non-rotational ly  symmetric 
perturbations.  His  calculations  showed  that  in  some 
cases  non-axi symmetri c modes  were  even  more  unstable  than 
the  axisymmetric  ones,  hence  it  was  important  to  take  them 
into  consideration.  He  also  found  that  the  surface  tension 
in  general  stabilized  the  fluid  motion,  except  for  axisym- 
metric perturbations  with  axial  wave  numbers  less  than  a 
certain  value.  This  result  was  first  shown  by  Rayleigh 
(1  879)  (also  Lamb  1932,  §273  ) . 

Chandrasekhar  (1960a,  b ; 1961,  §67,  §78b)  tried 
to  generalize  the  Ray  1 ei gh-Synge  criterion  for  arbitrary 
disturbances  and  for  rotating  flows  in  the  presence  of  an 
axial  flow  gradient.  He  came  to  the  conclusion  that  the 
Rayl ei gh-Synge  criterion  for  stability  continued  to  be 
valid  irrespective  of  perturbation  and  axial  flow  condi- 
tions. This  would  mean  that  a behavior  resembling  the 
Squire's  (1933)  transformation  for  plane  parallel  flows 
with  respect  to  three-dimensional  perturbations  could  be 
expected  in  the  swirling  flow  case.  However,  his  results 
appear  rather  implausible  because  of  the  characteristic 
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difference  between  the  axi symmetri c perturbati ons  for 
rotating  flows  and  the  non-axi  symmetri  c ones  for  swirling 
flows.  In  a mathematical  sense,  the  former  is  governed 
by  a Sturm-Liouvi  1 1 e system  where  the  principle  of  ex- 
change stabilities  is  usually  expected,  while  the  latter 
in  general  is  not.  In  a physical  sense,  the  former  con- 
cerns an  instability  of  the  Rayl ei gh -Tay 1 or  type,  in  which 
particle  exchange  is  restricted  in  axial  direction  (m  = 0) 
while  the  latter  is  an  instability  of  the  Kel vi n-Hel mhol tz 
together  with  Rayl ei gh-Tayl or  type,  in  which  particle  ex- 
changes occur  generally  along  the  direction  of  the  mean 
swirling  flows.  Howard  (1962)  first  pointed  out  that 
Chandresakhar ' s extension  of  Rayleigh's  criterion  was 
questionable  and  later  Fung  and  Kurzweg  (1973)  gave  an 
example  which  clearly  showed  that  Chandrasekhar's  result 
was  incorrect. 

Alterman  (1961)  studied  the  instabilities  of  cy- 
lindrical liquid  jets  bounded  by  a surrounding  liquid  and 
subject  to  rotation.  He  found  that  the  rotation-induced 
instability  for  axisymmetric  disturbances  depended  on  the 
angular  velocities  of  the  jet  relative  to  its  surroundings 

Batchelor  and  Gill  (1962)  made  a general  mathe- 
matical analysis  of  the  stability  of  axisymmetric  jets 
for  homogeneous  unbounded  fluids.  A condition  similar  to 
Rayleigh's  (1880)  inflection  point  theory  was  shown  to  be 
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necessary  for  instability  and  also  a semi-circle  theorem 
was  found  to  offer  a bound  for  possible  unstable  waves. 

Following  an  important  paper  by  Miles  (1961)  and 
a note  by  Howard  (1961)  on  the  stability  of  heterogeneous 
shear  flows,  Howard  and  Gupta  (1962)  obtained  two  stabil- 
ity conditions  for  axisymmetric  disturbances;  one  for  pure 
axial  flow  with  radially  dependent  density  and  another  for 
swirling  flows  with  constant  density.  In  both  cases  they 
found  that  a sufficient  condition  for  stability  required 
that  a non-dimensional  quantity,  reminiscent  of  the  Richard- 
son number,  be  greater  than  1/4  anywhere  within  the  flow 
region,  and  furthermore,  that  a semi-circle  theorem  held 
for  "statically  stable"  profiles  where  the  Rayl ei gh-Synge 
criterion  was  satisfied.  They  also  considered  flow  stabil- 
ity against  non-axisymmetric  disturbances,  but  no  general 
conclusion  could  be  drawn  from  their  stability  equation  in 
this  case  as  the  sufficiency  condition  was  found  to  be  vio- 
lated for  sufficiently  small  axial  wave  numbers.  This  re- 
sult does,  however,  suggest  that  azimuthally  periodic  per- 
turbations are  more  unstable. 

In  a paper  concerning  two-dimensional  vQrtex-type 
flows  with  constant  density,  Michalke  and  Timme  (1967)  de- 
rived a stability  equation  in  terms  of  perturbation  vortic- 
i ty . Two  analytical  and  one  numerical  solutions  were  found 
for  three  types  of  vortex  flows.  Their  results  revealed 
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that  the  two-dimensional  azimuthally  periodic  disturbances 
were  apparently  amplified  more  strongly  than  the  corres- 
ponding three-dimensional  ones.  However,  in  their  first 
example  of  the  cylindrical  vortex  sheet,  they  failed  to 
consider  the  centrifugal  effects  due  to  the  discontinuity 
of  velocity  at  the  interface,  and  accordingly  their  result 
must  be  modified.  An  analytical  example  by  Fung  and  Kurz- 
weg  (1973)  also  showed  that  the  flow  was  less  stable  with 
respect  to  azimuthally  periodic  disturbances  than  to  the 
axisymmetric  ones. 

Weske  and  Rankin  (1963)  investigated  experimentally 
the  motion  in  the  region  of  the  core  of  a vortex  which  ex- 
hibited peripheral  vorticity.  The  early  stages  of  instabil- 
ity observed  are  essentially  two-dimensional  with  large 
azimuthal  wave  numbers. 

Recently,  Johnston ( 1 972 ) observed  one  type  of  in- 
stability in  rotating  flows  by  injecting  a light  gas  tan- 
gentially into  a stationary  gas  of  higher  density.  The 
onset  of  instability  is  characterized  by  the  two-dimensional 
azimuthally  periodic  modes  which  are  apparently  most  unstable 
(see  Kurzweg  1 974) . 

Leibovich  (1969)  generalized  Howard  and  Gupta's 
(1962)  stability  condition  to  swirling  flows  with  radially 
dependent  density  but  restricted  his  discussions  to  axisym- 
metric perturbations. 
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Kurzweg  (1969)  derived  a general  stability  equa- 
tion for  heterogeneous  swirling  flows  with  radially  depen- 
dent velocity  and  density  for  arbitrary  infinitesimal 
disturbances  in  the  absence  of  viscous  and  gravitational 
effects.  A sufficient  condition  for  stability  was  found 
and  in  particular  it  was  shown  that  such  heterogeneous 
swirling  flows  could  be  guaranteed  stable  against  arbitrary 
infinitesimal  disturbances  when  the  density  was  an  increas- 
ing function  of  radius  and  at  the  same  time  the  velocity 
gradients  remained  small.  Later  Kurzweg  (1970)  obtained 
the  bounds  on  growth  rate  and  phase  velocity  for  the  in- 
viscid  ability  of  a constant  density  swirling  flow  with 
the  narrow  gap  approximation.  Two  types  of  swirling  pro- 
file were  investigated  numerically  and  the  effect  of  add- 
ing an  axial  velocity  to  the  Couette  flow  was  also  discussed. 

Here,  we  would  like  to  generally  investigate  the 
stability  characteristics  of  heterogeneous  swirling  flows 
and  compare  our  results  with  the  existing  theoretical  and 
ex  peri  mental  ones . 


CHAPTER  2 


MATHEMATICAL  FORMULATION 

We  wish  to  examine  the  stability  of  swirling 
flows  between  two  concentric  cyl  i nders  [R-j  , R2]  having 
the  steady-state  radially  dependent  velocity  field 
V [0,  rf2(r),  W (r)].  The  fluid  under  consideration  is 
assumed  inviscid,  incompressible  and  non-heat-conducting 
but  has  a radially  dependent  density  p^(r).  Euler's  equa- 
tions of  motion,  the  incompressibility  condition  and  the 
continuity  equation  governing  such  a swirling  flow  in 
cylindrical  co-ordinates  (r,  6,  z)  are,  respectively, 

(||  + V • VV)  = - -VP  + g (2.1) 

O L P 

(2.2) 

(2.3) 

To  examine  the  stability  of  this  heterogeneous  flow, 
we  postulate  an  initial  set  of  solutions  based  on  the  summa- 
tion of  the  steady-state  flow  and  infinitesimal  time-depen- 
dent  fluctuations,  namely. 


IP 

3t 


^ + V 


(Vp) 


V = 0 
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= 0(r,e,z;t) 

Vq  = rn(r)  + ^(r,0,z;t) 
= W(r)  + «(r,9,z;t) 

P = P^{r)  + p{r,0,z;t) 

p = pQ^*")  ^(>">0.2;^) 


(2.4) 


Substituting  equations  (2.4)  into  equations  (2.1),  (2.2) 
and  (2.3),  neglecting  the  gravitational  effects,  and  fol- 
lowing the  normal  linearization  procedure  of  dropping 
second  and  higher  order  terms  in  the  assumed  infinitesimal 
perturbations  lead  to  the  linearized  differential  equations 


+ PQ(d*u  - 2Q9)  - = 0 

II  + Pgr[d*9  + D*(rf^)Q]  = 0 

If  + Po[d*(i  + (DW)Q]  = 0 
d*p  + (Dp^)Q  = 0 


M + 0 + 1 ii 

9 r r r 90 


9z 


L 

(2.5a) 

(2.5b) 


( 

(2.5c) 


(2.5d) 


(2.5e) 
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where 


D 


dr  , 


* 1 
D = D + - 
r 


and 


A +qA. 

9t  90 


+ w 


9z 


with  boundary  conditions  d = 0 at  R-|  and  R2. 


Eliminating  Q,  W and  p from  equations  (2.5),  we  obtain 
the  following  two  equivalent  partial  differential  equa- 
tions for  stability. 


1 

2 

r 


9^p  9^ 

90^  9z^ 


Po[d*(D*a) 


P A _ (DW)  ii] 

r 90  9z-' 


0 (2.6) 


r 90 


+ PQ^*P 


+ 4>u  = 0 
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C 

(2.7) 


where 


3 D[Pn(r^^)^] 

r3 
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The  coefficients  of  equations  (2.5)  or  equations  (2.6) 
and  (2.7)  depend  only  on  the  radial  coordinate  r,  while 
the  remaining  independent  variables  z,  6 and  t enter  only 
in  the  partial  derivatives.  This  allows  one  to  assume 
that  (2.5),  (2.6)  and  (2.7)  admit  solutions  having  an 
exponential  behavior  in  z,  9 , and  t.  This  assumption  re- 
duces all  these  partial  differential  equations  to  ordinary 
differential  equations. 

Following  the  normal  mode  analysis  indicated  in 
Chapter  1,  we  apply  the  Laplace  transform  with  respect  to 
time  t and  the  Fourier  transform  with  respect  to  the  spa- 
tial variables  (z,  9).  The  transforms  have  the  form 

{p,u}  (k,m,0);r)  = f f f { p ,u  } ( r , z , 9 ; t ) e^  ^"^‘™'“^^dzd9dt 

0 0 

(2.8) 


which  is  equivalent  to  the  standard  transformation  of  assum- 
i ng 


{P.Q} 


{p(r) ,u(r)}e 


i ( kz+m9 -u)t ) 


We  are  here  considering  only  temporal  disturbances  so  that 
the  axial  wave  number  k is  real  while  the  amplification 
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factor  is  complex;  m represents  the  azimuthal 

wave  number  and  takes  on  only  integer  values  because  of 
the  periodicity  of  the  solution  in  the  6 direction.  Using 
transform  (2.8)  in  equations  (2.6)  and  (2.7)  leads  to  the 
two  ordinary  differential  equations 


(k^  + = ipQ{[k(DW)  + niO-^^]u  - N(D*u)}  (2.9) 

r ^ 

NDp  + ^p  = i($  - PqN^)u  (2.10) 


By  eliminating  p,  these  equations  may  be  reduced,  after 
some  manipulation,  to  the  single  second  order  stability 
equa  ti on 


N^[D(PqED*)  - pQk^]u  + u[NF  + G]  = 0 


(2.11) 


wi 


ith  boundary  conditions  li(R.|)  = u(R2)  = 0 


Here 


N = kW  + mSi!  - (1) 


s = 


m 

kr 
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E = (1  + s^)-'' 

F = -kD*{p^E[DW  + sD*(rf2)]} 

G = k^E[1>  - 2sPqJ^(DW)  + s^rf2^(Dp^)] 

D*=  D - - 

Equation  (2.11)  was  first  derived  by  Kurzweg  (1969) 
and  reduces  to  the  one  given  by  Howard  and  Gupta  in  the 
limiting  case  of  constant  density.  Together  with  the 
homogeneous  boundary  conditions,  equation  (2.11)  forms  an 
eigenvalue  problem  for  the  complex  eigenvalue  w for  given 
specific  values  of  m,  k,  W and  p^ . The  system  is  un- 
stable whenever  > 0.  However,  one  should  keep  in  mind 

that  since  (2.11)  is  invariant  under  complex  conjugation, 
any  solution  with  < 0 (damped  wave)  implies  the  exis- 
tence of  another  solution  with  > 0 (growing  wave),  and 
hence  instability.  This  fact  shows  that  in  the  absence 
of  viscous  effects  there  can  be  no  damped  wave  and  the  flow 
is  either  unstable  (w^.  > 0)  or  neutrally  stable  (oo^-  = 0). 

In  order  to  distinguish  the  bounding  neutral  curve  (w^.  = 0) 
characterizing  stability  from  the  other  neutral  curves 
(also  0)^.  = 0)  for  all  other  stable  solutions,  we  define, 
following  Miles  (1961),  the  neutral  curve  as  a stability 
boundary  if  and  only  if  there  exist  contiguous  eigenvalues 
for  which  > 0 (unstable  modes). 


CHAPTER  3 


SOME  GENERAL  STABILITY  CRITERIA 

After  Rayleigh's  (1916)  derivation  of  his  condi- 
tion for  stability  of  rotating  flows  subject  to  axisym- 
metric  disturbances  using  a conservation  of  angular 
momentum  argument,  Synge  (1933)  proved  the  same  stability 
condition  for  heterogeneous  fluids  by  casting  the  problem 
into  the  form  of  a Sturm-Liouvi  1 1 e system  and  using  the 
properties  of  this  system.  The  resultant  well-known 
Rayleigh-Synge  criterion  obtained  is  that  the  flow  can  be 
guaranteed  stable  for  infinitesimal  axisymmetric  distur- 
bances whenever  the  product  of  its  density  and  the  square 
of  its  circulation  is  an  increasing  function  of  radius. 

We  here  give  an  alternate  derivation  of  the  Ray- 
leigh- Synge  criterion  following  Synge's  approach  but  without 
assuming  that  the  perturbations  are  periodic  in  the  z direc- 
tion. 

For  symmetry  about  the  axis  and  zero  axial  flow, 
equations  (2.5)  can  be  reduced  to 

II  ^o^It  ■ = 0 (3.1a) 
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9t 


+ D*(rJ^)(i  = 0 


(3.1b) 


az 


(3.1c) 


(DPq)u  = 0 


(3. Id) 


M + i + M 

9r  r 9z 


0 


(3.1e) 


Equations  (3.1a)  through  (3. Id)  can  be  combined  into  the 
following  partial  differential  equation 


(3.2) 


A. 

A1  SO  we  define  a Stoke's  stream  function  'F  such  that 


0 


i ii 

r 9z 


i 

r 9 r 


(3.3) 


the  continuity  equation  (3.1e)  will  be  satisfied  automati- 
cally, and  (3.2)  can  be  transformed  to 
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2 2^  2 D ^ 

+ X)  ill  + ^ {JL  il)}=  0 

‘n2  •’o  9Z=^  ‘p„ 


(3.4) 


Letting  the  infinitesimal  perturbations  vary  exponentially 
in  time  (Yih  1965,  Chapter  6)  such  that 


¥ = Y{r  ,z)e 


-ioat 


allows  us  to  reduce  (3.4)  to  the  form 


9^'i'  + n ii  = 0 

9r  CO  0 9z  0 


(3.5) 


34' 


with  boundary  conditions  ^ ^ at  R-j  and  R2 


The  equation  for  characteristics  of  (3.5)  is 


^ = ± /4(-^)-  1 
‘‘P  Po 


(3.6) 


1 , j, 


Equation  3.5  is  hyperbolic  if  -2^”^“^  (3.7a) 


(A)  ^0 


0.2  Pq 


is  parabolic  if 


(3.7b) 
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and  is  elliptic  if  -y(— ) -i  <0  (3.7c) 

There  are  two  families  of  real  characteristic 
curves  describing  the  physical  characteristic  surfaces 
in  the  hyperbolic  case,  only  one  in  the  parabolic  case 
and  none  in  the  elliptic  case.  For  this  reason,  two 
conclusions  can  be  drawn  just  by  discussing  the  charac- 
teristics of  the  partial  differential  equation  (3.5). 

O 

First,  real  characteristics  do  not  exist  if  w 

is  complex  since  is  always  real.  In  order  to  have  at 

^0 

least  one  set  of  real  characteristic  curves  satisfy  a 

O 

physical  situation,  it  requires  co  be  real,  i.e.,  the 
amplification  factor  is  either  real  (stability)  or  pure 
imaginary  (instability).  Hence  the  principle  of  exchange 
stabilities  holds  for  the  flow  under  consideration  with 
to  = 0 corresponding  to  the  neutral  stability  condition. 

Second,  the  assumed  stream  function  T (r,  z)  is  a 
scalar  function  such  that 


d^ 


9r 


dr  + 


3z 


dz 


0 


9r 


or 


dz 
d r 


3'F 

Dz 


(3.8) 


22 


Substituting  (3.3)  for  exponential  solutions  in  time,  i.e., 

{Q,  w)  = {u(r,z),  w(r,z)} 

in  (3.8)  yields 


w 


Figure  3.1  Fluid  particle  movements  in  the 
velocity  field  for  ax i symmetric 
di sturbances. 

From  the  definition  of  the  stream  function,  the  fluid 
particles  should  move  along  the  streamlines  described  by 
the  stream  function  with  the  slope  governed  by  equation 
(3.9)  (see  Figure  3.1).  In  the  case  under  consideration 
here  the  equation  describing  the  slope  of  the  stream  func- 
tion happens  to  be  the  equation  for  characteristics  of  the 
partial  differential  equation  (3.5)  governing  stability. 

In  order  to  have  real  streamlines  as  well  as  real  fluctuation 
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velocities  u and  w exist,  it  is  required  that  the  slope 
of  streamlines  in  (3.6)  and  (3.9)  be  real,  i.e.,  equation 
(3.5)  can  either  be  hyperbolic  or  parabolic.  Accordingly 
the  criterion  governing  stability  is 

-l(— )-l  > 0 (3-10) 

J Po 

2 . . 2 

i.e.,  $ > PqW  for  positive  w , 

stabi 1 i ty 

2 . 2 

and  $ < PqW  for  negative  w , 

i n stabi 1 i ty 


The  well-known  Rayl eigh-Synge  criterion,  $>o,  then  fol- 
lows from  the  principle  of  exchange  stabilities  with<^=  0 
corresponding  to  neutral  stability. 

Perhaps  we  should  remember  that  this  is  the  criter- 
\ 

ion  for  instability  of  the  Rayl eigh-Tayl or  type  since  the 
mechanism  is  based  on  the  particle  exchanging  due  to  the 
centrifugal  (or  circulation)  effects.  It  should  be  empha- 
sized once  more  that  this  criterion  applies  only  to  axisym- 
metric  perturbations  with  zero  axial  velocity  (or  constant 
axial  flows).  Once  perturbations  are  allowed  to  be  non- 
axisymmetric  or  the  axial  velocity  gradients  are  taken  into 
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consideration,  new  mechanisms  for  instability  will  arise 
and  instabilities  of  the  Kel vi n-Helmhol tz  type  can  then 
take  place. 

The  investigations  concerning  instabilities  of 
the  Kel vi n-Helmhol tz  type  for  heterogeneous  shear  flows 
was  first  attempted  by  Taylor  (1931)  and  Goldstein  (1931). 

By  discussing  the  complex  zeros  of  the  Bessel  functions, 
they  discovered  that  the  Richardson  number,  defined  as  the 
ratio  of  the  density  gradient  of  the  flow  to  its  shearing 
rate,  equal  to  1/4  was  a critical  value  for  stability. 

(The  governing  stability  equation  was  then  proposed  by 
Thorpe  (1969)  to  be  called  the  Tayl or-Gol dstei n equation, 
in  recognition  of  its  first  use  by  Taylor  and  Goldstein 
in  determining  the  stability  of  certain  stratified  flows.) 
Later  Miles  (1961)  considered  a more  general  stratified 
flow  and  proved  the  same  sufficiency  condition  by  applying 
the  method  of  Frobenius  to  the  Tayl or-Gol dstei n equation, 
with  a restriction  that  the  velocity  gradient  was  nowhere 
equal  to  zero,  which  was  removed  by  Howard  (1961).  In  addi- 
tion, Howard  illuminated  a way  to  the  study  of  shear  in- 
stabilities by  proving  the  well-known  semi-circle  theorem, 
which  has  been  broadly  used  in  the  theoretical  investiga- 
tions of  parallel  flows  as  well  as  swirling  flows. 

Howard  and  Gupta  (1962)  published  an  important 
paper  concerning  stability  of  swirling  flows  without  density 
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variations,  in  which  instabilities  of  the  Kel vi n-Hel mhol tz 
type  were  discussed.  We  now  want  to  examine  the  general 
features  of  stability  for  heterogeneous  swirling  flows. 

Two  types  of  instabilities  are  under  consideration:  "cen- 

trifugal instabilities"  representing  unstable  situations 
caused  by  centrifugal  forces;  and  "shear  instabilities" 
corresponding  to  unstable  solutions  produced  by  shearing 
effects.  The  former  is  an  instability  in  the  Rayleigh- 
Taylor  sense  while  the  latter  belongs  to  an  instability 
of  the  Kel vi n-Hel mhol tz  type. 

Following  Howard  (1961),  Howard  and  Gupta  (1962), 
we  assume  there  exists  an  unstable  mode  (w^.  > 0)  and  let 
u = where  ? is  a real  constant.  Substituting  this 

in  equation  (2.11),  one  obtains 

D[(pQEN^^)D*ia  + + N[F  + (;D*(PqEDN)]  + 

[G  + c(c  - 1)PqE(DN)^]}i|^  = 0 (3.11) 

Multiplying  (3.11)  by  rijj  (the  complex  conjugate  ofr^J^) 
and  integrating  over  the  flow  region,  we  make  use  of  the 
boundary  conditions  iJ^(Ri)  = iKR2)  =0  to  obtain  the 
integral  form  of  the  stability  equation 

/ N^^"^{-P  + k^UI^]  + N[F  + cD*(p  EDN)]|(p|^  + 

R^  ^ ° 

[G+  (.{(.  - 1 )()^E(DN)^]  I 'H  ^ 1 rdr  = 0 


(3.12) 
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Several  stability  criteria  for  the  heterogeneous 
swirling  flows  under  consideration  can  be  derived  from 
this  integral  representation.  These  are  now  obtained  as 
fol 1 ows : 

3 . 1 Sufficiency  Condition  for  Stability 

Letting  c=  (3.12)  leads  to  the  integral 

«2 

/ {-PqN[E|D%|2  + + {[F  + D*(pqEDN)]  + 

•^1 

^-{G  - -J;  PoE(DN)^]}|/}  rdr  = 0 (3.13) 


whose  imaginary  part  is 

* I 1 2 

‘^i^  + [G  - I PQE(DN)^]j^}rdr  = 0 (3.14) 


It  can  be  seen  from  (3.14)  that  unstable  solutions  (oo  ^ > 0) 
are  not  possible  if  the  integrand  is  positive  definite 
throughout  the  whole  integration  domain.  Therefore  the 
flow  can  be  guaranteed  stable  if  the  function  { G-  ^PqE(DN  )} 
is  positive  everywhere  inside  the  flow  region,  that  is. 
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- ^-(Dwft  - 2s(DW)[Q  + + s^[rJ^^D(log  p^)  - ^ r^(Df^)^]  > 0 

(3.15) 


This  criterion  was  first  obtained  by  Kurzweg  (1969)  and 
is  reducible  to  Howard  and  Gupta's  (1962)  result  for  con- 
stant density. 

An  alternate  derivation  of  condition  (3.15)  can  also 
be  obtained  by  considering  a Frobenius  solution  to  equa- 
tion (2.11)  following  an  approach  used  by  Miles  (1961)  in 
a related  study  on  stratified  shear  flows  • The  derivation 
is  given  as  follows: 

Letting  u = ^ allows  one  to  transform  equation 
(2.11)  to  the  form 


D (f)  + 


D*(P  E) 


-Dcj)  + 


Po^‘ 


= 0 (3.16) 


The  critical  point  r - r^  of  this  equation  is  defined  by 
N = 0 

or  equivalently  kW  + mJ^  = co 


with  a restriction 


DW  + srDJ^ 
r 


in  (Ri , R^) 


f 0 


It  is  a regular  singularity  for  the  differential  equation, 
and  leads  to  an  indicial  equation 


- K + g = 0 
^ c 


or 

9'c  <3.17) 

where  o o 

$ - 2sp  Jl(DW)  + s rfl^(Dp^) 

g = P 

^ Pq(DW  + srDn)"^ 


The  two  linearly  independent  solutions  to  equation  (3.16) 
are  in  the  form 


(r) 


9c 


5 c,(r-r Jj  (3,18) 


where  C^.  is  arbitrary  constant.  The  exponent  of  the  solu- 
tions will  be  complex  whenever  g^  > ^,  a condition  identical 
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with  that  given  by  (3.15),  Perhaps,  one  should  go  through 
some  additional  mathematical  arguments  in  deriving  this 
result  by  considering  the  Reynolds  stress  as  done  by  Miles 
(1961).  However,  it  follows  from  the  sufficiency  condition 
proven  by  the  integral  technique  that  neither  (j)+  nor  <j)_ 
may  have  complex  zeros  and  still  match  the  boundary  condi- 
tions when  ^ • The  upper  bound  on  which  possible  un- 
stable solutions  can  ever  exist  is  thus  g^  = . 

We  have  just  shown  from  the  above  analysis  that  the 
sufficiency  condition  for  stability  depends  completely  on 
the  second  order  singularity  of  the  differential  equation 
(2.11).  As  a result  of  it,  two  phenomena  concerning  stabil- 
ity analysis  can  be  seen  immediately: 

First,  neither  the  neglect  of  the  inertia  effects  of 
density  variation  nor  the  narrow  gap  analysis  will  change 
the  second  order  singularity  of  the  stability  equation,  so 
that  the  sufficiency  condition  will  remain  unchanged.  (A 
similar  argument  can  be  applied  to  the  familiar  Taylor- 
Goldstein  equation  arising  in  the  stability  of  stratified 
shear  flows. ) 

Second,  when  Howard  and  Gupta  (1962)  made  this  in- 
vestigation of  homogeneous  swirling  flow  stability  they 
realized  that  their  sufficiency  condition  for  stability 
was  always  violated  for  sufficiently  small  axial  wave 
numbers  and  no  general  stability  criterion  was  obtainable. 
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In  order  to  explain  this  situation,  one  should  look 
back  at  the  original  differential  equation  governing 
stability.  For  constant  density  and  k = 0,  the  second 
order  singularity  of  equation  (2.11)  will  be  eliminated 
and  the  equation  they  derived  becomes  an  equation  of 
the  Rayleigh  type  having  no  second  order  singularity. 
Another  look  at  the  Frobenius  solution  (3.18)  shows 
that  the  exponents  of  the  indicial  equation  (3.17)  are 
always  real  and  equal  to  an  integer  because  of  the  ab- 
sence of  the  second  order  singularity.  The  second  lin- 
early independent  solution  is  then  proportional  to 
<j)^(r)log(r-r^)  and  the  solution  for  the  perturbation 
velocity  described  by  (3.18)  may  possess  a complex  zero 
matching  the  boundary  conditions.  Under  this  situation, 
the  possibility  of  non-trivial  solutions  for  unstable 
modes  always  exists  and  no  stabilities  can  then  be  guaran- 
teed before  solving  the  equation.  For  this  reason,  the 
sufficiency  condition  derived  by  Howard  and  Gupta  fails  to 
predict  stabilities  due  to  the  disappearance  of  the  second 
order  singularity  in  their  equation  for  sufficiently  small 
k.  (Based  on  the  same  argument,  the  Richardson  number 
criterion  in  the  Taylor-Goldstein  equation  becomes  mean- 
ingless in  Rayleigh's  equation  for  homogeneous  fluids 
since  the  Richardson  number  is  identically  equal  to  zero 
due  to  the  absence  of  buoyancy  effects,  and  no  general 
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information  can  be  obtained  as  a result  of  this  condi- 
tion. ) 

In  fact,  the  physical  meaning  of  the  sufficiency 
condition  which  is  in  terms  of  the  second  order  singular- 
ity for  the  equations  governing  stability  represents,  in 
some  way,  the  potential  stabilizing  effects  existing  in 
the  flow  to  suppress  instability.  This  phenomenon  can 
be  explained  as  follows: 

Consider  two  particles  A and  B with  equal  fluid 
volumes  in  a two-dimensional  parallel  velocity  field 
shown  in  Figure  3.2. 


Assume  a statically  stable  profile  with  < pg  in  a 
heterogeneous  density  field.  The  interchange  of  position 
between  particle  A and  B is  possible  only  when  the 


^ U 


Figure  3.2  Fluid  particle  interchange  in 
the  gravitational  force  field. 
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interchange  kinetic  energy  gained  by  these  two  particles 
is  sufficiently  large  to  overcome  the  work  done  due  to  the 
exchange  of  them  in  the  gravitational  field.  As  a result 
of  it,  stability  can  be  assured  whenever  the  required  work 
done  exceeds  the  change  in  kinetic  energy.  In  other  words, 
instabilities  are  not  possible  when  stabilizing  effects  are 
strong  enough  to  suppress  instabilities  due  to  shear. 

However,  the  interchange  of  the  same  two  particles 
in  a velocity  field  without  gravitational  effects  or  in 
a homogeneous  density  field  with  gravitational  effects 
will  result  in  no  work  done  or  potential  energy  changed. 

Due  to  the  absence  of  stabilizing  effects,  any  change  in 
kinetic  energy  may  cause  instability.  This  is  why  the 
Richardson  number  criterion  fails  to  predict  stability 
for  homogeneous  fluids. 

The  same  arguments  can  be  used  to  explain  the 
stability  for  non-rotating  jets  and  for  the  homogeneous 
swirling  flows  against  azimuthally  periodic  perturbations, 
with  the  centrifugal  force  field  replacing  the  gravita- 
tional field.  In  both  cases,  no  stabilizing  factors  exist 
and  any  axial  or  rotational  shear  component  may  upset  the 
stability  of  the  flow,  so  that  no  general  conclusion  for 
stability  or  possible  unstable  waves  can  be  drawn.  In  the 
latter  case,  it  suggests  that  the  flow  is  less  stable  to 
two-dimensional  azimuthally  periodic  modes  (i.e.,  k = 0) 
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than  the  corresponding  three-dimensional  ones  with  non- 
zero axial  wave  numbers. 

Condition  (3.15)  can  be  expressed  in  terms  of  a 
non-dimensional  local  quantity  defined  as 

X = ^ (3.19a) 

p^[DW  + s(rDfi  + 4f2)] 

If  DW  + s(rDil  + 4^1)  = 0 somewhere,  we  allow  A ^ . In 

terms  of  this  quotient  the  stability  condition  given  by 
(3.15)  is 

X J I (3.19b) 

for  all  r in  (R^,  R2).  The  result  obtained  by  Leibovich 
is  deducible  from  this  equation  for  the  special  case  of 
ax i symmetric  perturbations  (i.e.,  s = 0).  Also  the  well- 
known  Rayleigh-Synge  criterion  can  be  regained  from  (3.19b) 
in  the  absence  of  axial  flow  gradients  and  for  s = 0. 

This  non-dimensional  number  A,  reminiscent  of  the 
Richardson  number  criterion  for  stability  of  heterogeneous 
shear  flows,  is  defined  as  the  ratio  of  centrifugal  effects 
to  shearing  effects.  Condition  (3.19b)  can  be  thought  of 
as  the  criterion  for  shear  stability  in  flows  with  "cen- 
trifugally  stable"  profiles.  Here  we  say  that  a flow  is 
"centrifugally  stable"  if  the  product  of  its  density  and 
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the  square  of  its  circulation  is  nowhere  a decreasing 
function  of  radius.  Condition  (3.19b)  clearly  reveals 
that  the  centrifugal  stability  (i.e.  $>0)  does  not  neces- 
sarily lead  to  stability  for  nonaxi symmetri c perturbations 
or  in  the  presence  of  axial  flow  gradients.  Accordingly 
it  is  concluded  that  the  Rayleigh- Synge  criterion,  a condi- 
tion for  centrifugal  stability,  is  neither  necessary  nor 
sufficient  for  shear  stability. 

It  should  be  noted  that  Chandrasekhar  in  his  book 
fails  to  realize  this  point  and  proposes  to  show  that  <!>>0 
is  a necessary  and  sufficient  condition  for  the  stability 
of  swirling  flows  subject  to  arbitrary  perturbations.  As 
pointed  out  by  Howard  (1962)  such  a criterion  is  certainly 
incorrect. 

Violations  of  the  Rayleigh- Synge  condition  imply 
centrifugal  instabilities  for  axisymmetric  perturbations 
and,  at  the  same  time,  result  in  the  violation  of  the  suf- 
ficiency condition  (3.19b),  which  may  induce  shear  instabil- 
ities (even  though  it  is  not  necessarily  so).  Hence,  two 
types  of  instabilities  may  be  expected  to  take  place  simul- 
taneously whenever  <!><0. 

Here  we  wish  to  point  out  that  equation  (2.11)  forthe 
narrow  gap  approximation,  i.e.,  R^  - ■*'  ^2^’  ^ 

the  absence  of  axial  velocity  gradients  ( DW  = 0)  is  mathe- 
matically similar  to  the  Tayl or-Goldstei n equation  encountered 
in  two-dimensional  stratified  shear  flows. 
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Setting  ^ = 1 allows  us  to  rewrite  equation  (3.11) 
as 

D[(PqEN^)D%]  + {-ppkV  + N[F  + D*(PqEDN)]  + G}(*;  = 0 (3.20) 

Letting  r = Rq  + y and  considering  the  small  gap  approxi- 
mation y « R 
^ 0 


so  that 


wi  th 


* ri 


s s 


m 

0 " kR. 


and 


*^o  " 2^^1  ^ ^2^ 


allow  us  to  reduce  (3.20),  after  neglecting  the  axial 
velocity  gradients,  to  a relatively  simple  form 


^[PqOJ  - e)^^]  + [G  - P„k(l  + s^)2(U  - 8)]*  = 0 (3.21) 


wi  th 


U(y)  = rSl 


?„(f  - “o> 


1 


H i|2  dp  ,.2 

d/_U\  , 0/U^ 

2 dy^PoR„^  dy  ^R„^ 
S 0 0 
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If  the  centrifugal  term  produced  by  the  steady-state 
rotation,  has  only  small  variations  across  the  gap,  then 


G 


1 ii2  dp 

( 1 + -L)  (— ) — - 

S O'' 

0 


and  equation  (3.21)  takes  the  form 


.2  dp, 


^[p„(U  - c)2g]  + [(1  + - p/(1  ^ s2)(U  - c')2]*  = 0 (3.22) 


dy 


This  equation  is  mathematically  identical  to  the  Taylor- 

Goldstein  equation  before  the  effects  of  density  variation 

in  the  inertia  term  are  neglected,  provided  the  centrifugal 

force  termCn— } is  replaced  by  the  gravitational  constant 
^0 

g.  The  sufficiency  condition  for  stability  of  rotating 
flow  in  this  narrow  gap  approximation  is  found  to  be 

(1  + -L)J  = [1  + (-;r)  ^ I (3.23) 

wi  th 

d 


dy 

^0  dy ^ 


- p ' g 
o^ 

P.(U 


being  the  well-known  Richardson  number  encountered  in 
heterogeneous  parallel  flows.  Equation  (3.23)  will  have 
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a minimum  for  each  local  "modified  Richardson  number" 

[1  + (kR^/m)  ]J  or  for  each  given  velocity  and  density  pro- 
file when  k goes  to  zero,  or  m goes  to  infinity,  implying 
that  the  flow  is  less  stable  with  respect  to  two-dimensional 
azimuthally  periodic  perturbations  than  to  the  corresponding 
three-dimensional  ones  with  k 0.  This  in  turn  would  imply 
that  the  Squire  (1933)  theorem  (see  also  Yih  1955)  for  plane 
parallel  flows  with  respect  to  three-dimensional  perturba- 
tions would  be  expected  in  revolving  flows  as  well,  at  least 
for  large  azimuthal  wave  numbers.  This  argument  is  supported 
by  the  analytical  results  by  Michalke  and  Timme  (1967)  and 
by  those  of  Fung  and  Kurzweg  (1973),  as  well  as  experimental 
results  by  Weske  and  Rankin  (1963)  and  by  Johnson  (1972)  (see 
also  Kurzweg  1974). 

Following  Taylor  (1931)  and  Goldstein  (1931),  Drazin 
(1958)  showed  analytically  that  a quarter  was  the  minimum 
value  of  the  Richardson  number  J for  which  the  flow  can  be 
guaranteed  stable  for  stratified  shear  profiles.  Later  on 
Hazel  (1972)  confirmed  this  point  by  numerical  investiga- 
tions of  the  Taylor-Goldstein  equation.  Perhaps  one  may 
question  if  the  sufficiency  condition  ^ ^ for  hetero- 

geneous swirling  flows  can  be  improved.  That  this  gener- 
ally cannot  be  done  is  demonstrable  by  considering  a 
homogeneous  flow  for  rigid  body  rotation  subject  to  azi- 
muthally periodic  disturbances.  This  flow  is  neutrally 
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stable  (as  will  be  shown  through  detailed  calculation 
later)  at  precisely  X = ^. 

3 . 2 Necessary  Condition  for  Instability 

One  of  the  necessary  conditions  for  centrifugal 
instability  subject  to  ax i symmetric  disturbances  is  well 
described  by  the  Rayleigh- Synge  criterion  for  zero  axial 
velocity  gradients.  Two  other  necessary  conditions  for 
shear  instability  can  be  derived  from  the  general  integral 
equation  (3.12)  by  letting  ; = 0.  This  yields 

/ {-PqCEID  i>\^  + + [^  + -^]  = 0 (3.24a) 

R-|  N 

whose  real  part  is 

*^2  2 
f {F|N|^(kW  + mf^  - to^)  + G[(kW  + mQ  - w^)^  - w^]}|-!^^rdr  = 

/Pq[E|D*4^|^  + k^|i|;|^]rdr  (3.24b) 

R T 


and  imaginary  part  is 
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0).  / {F  |N  1^  + 2G(kW  + tnfi  - w 

R, 


(3.24c) 


We  have  not  succeeded  in  drawing  any  general  conclusion 
from  (3.24c)  for  possible  unstable  solutions;  however, 
some  special  cases  in  which  the  function  (G)  vanishes  can 
be  obtained.  In  such  a case,  exactly  as  in  the  proof  of 
Rayleigh's  (1880)  inflection  point  theorem  for  instability 
of  inviscid  parallel  flows,  we  can  conclude  that  the  func- 
tion inside  the  curly  bracket  in  (3.24c)  should  change 
sign  in  (R-|,  R2)-  Certain  special  cases  of  equations  (3.24) 
can  lead  to  the  following  results: 

3.2.1  Non-axi symmetri c Perturbations  of  Heterogeneous 


Axial  Flows 


Equations  (3.24b,c)  for  = 0 yield 


R-i  I W - C I K 


< 0 


(3.25a) 


c.  / D*(p^EDW) 
R 1 


rdr  = 0 


(3.25b) 
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where 


In  order  to  have  unstable  solutions  (c^  > 0)  we  require  that 

the  quantity  { _}  have  at  least  one  numerical  ex- 

m^  + ( kr ) 

tremum  value  at  some  point  in  (R-j,  R^).  This  is  exactly 
the  condition  Batchelor  and  Gill  (1962)  derived  except 
that  they  only  studied  the  homogeneous  jets  (see  also  Ray- 
leigh 1892a)L  However,  the  density  variation  in  no  way  in- 
fluences the  characteristics  of  the  general  criteria  (see 
both  the  necessary  condition  and  the  later  to  be  discussed 
semi-circle  theorem)  for  the  pure  axial  flows  if  the  gravi- 
tational effects  are  neglected. 

An  extension  of  Rayleigh's  theorem  by  Fjdrtoft 
(1950)  and  H0iland  (1953)  states  that  the  necessary  condi- 
tion for  instability  is  that  the  vorticity  should  somewhere 
have  a maximum  inside  the  flow  domain.  A similar  refine- 
ment of  the  above  necessary  condition  for  instability  of 
radially  dependent  jet  flows  can  be  also  made.  Let  r^  be 
the  point  where  the  function  {D*(p^EDW)}  vanishes,  the 
corresponding  axial  velocity  being  W^.  Equations  (3.25b) 
can  be  written  as 

R 2 2 

(W^  - c J / D*(p  EDW)— ^rdr  = 0 

^ R^  ° |W  - cr 


(3.25c) 
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Equations  (3,25a,c)  then  yield 


f (VI  - W )D*(p  EDW) — ^rdr  =-/p  [E|D  i>\^  + | | ^ ]rdr  < 0 

^ 0 |w  - cr  Ri 


From  this  follows  a stronger  necessary  condition  for  instabil- 


ity, namely  that  (W  - W^)D*(p^EDW)  < 0 somewhere  in  the  flow 
field.  In  particular,  if  W(r)  is  a monotone  function  and 
D*(p^EDW)  vanishes  only  once  in  (R^,  R2).  a necessary  condi- 
tion for  instability  is  that  (W  - W )D(  '))  ^ 0 

s'  + k^r^ 

through  the  flow  region. 


3.2.2  Azimuthally  Periodic  Perturbations  of 
Homogeneous  Swirling  Flows 

For  k = 0 and  Dp^  = 0,  the  sufficiency  condition 

(3.19b)  fails  to  predict  stability  because  of  the  absence 

of  centrifugal  stabilizing  effects.  The  necessary  condi- 

* 

tion  for  instability  is  that  DD  (rfi)  should  change  sign, 
i.e.,  the  vorticity  in  axial  direction  should  have  a maxi- 
mum in  (R.|,  R^).  This  result  is  due  to  Rayleigh  ( 1880) 
for  zero  axial  velocity;  the  presence  of  axial  flows  of 
course  has  no  relevant  effects  when  the  perturbations  are 
azimuthally  periodic.  Unquestionably  this  necessary  condi- 
tion can  also  be  refined  using  Fj^rtoft  and  Hfiiland's  method. 
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For  instability,  the  term  {(fl  - fl^)D*(rfi)}  should  be 
negative  somewhere  in  the  field  of  flow,  being  the 
value  of  the  angular  velocity  where  the  axial  vorticity 
component  vanishes. 

3 . 3 Bound  for  the  Complex  Wave  Speed 

A general  bound  on  the  growth  rate  for  centrifugal 
and/or  shear  instability  can  also  derived  from  the  general 
integral  equation  (3.12).  For  C = 1,  equation  (3.12)  has 
the  complex  form 

/ {N^p  [ElD*ipl^  + - {N[F  + D*(p^EDN)]  + G}  1 [ ^}rdr  = 0 

R.J  ^ ° 

(3.26) 


with  the  real  part 


/ { [ ( kW  + mS^  - 


- {(kW  + ma  - U|.)[F  + D*(p^EDN)]  + G)|*|^r)dr  = 0 (3.27a) 


and  imaginary  part 
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^2  ^2  ^2 

/(kW  + mfl)Qdr  = o)  /Qdr  + |-/[F  + D*  ( p^EDN  ) ] | ij;  | ^rdr  (3.27b) 
where 

Q = Pq[E|D  \l)\^  + k^|(|;|^]r  > 0 
Substituting  (3.27b)  to  (3.27a),  we  obtain 
R 2 R 2 R 2 

/ (kW  + mS^)^Qdr  = (o3^  + lo?)  / Qdr  + / {(kW  + mfi) 

R,  «1  «1 

•[F  + 0,(Pj|EDN)]  + G)  |*|^rdr  (3.27c) 


The  above  integrals  are  reminiscent  of  those  derived  by 
Howard  (1961),  and  Howard  and  Gupta  (1962)  in  their  formu- 
lation of  the  semi-circle  theorem  serving  as  the  bound  of 
the  complex  velocity  in  stratified  shear  flows. 

Two  types  of  semi-circle  bound  for  instability  of 
swirling  flows  will  now  be  discussed. 

3.3.1  The  Semi-circle  Theorem  for  Shear  Instabilities 

For  the  present  swirling  flows,  the  semi-circle 
theorem  can  also  be  derived  for  special  profiles.  We 
consider  the  following  four  categories; 
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(a)  For  non-rotating  jet  flows,  = 0 , equations 
(3.27b,c)  become 


R 


2 

/ WQdr 


R 


1 


/ Qdr 


(3.28a) 


/ W^Qdr  = 


R 


1 


(c^  + c^)  / Qdr 
R] 


(3.28b) 


where  c = ^,  the  complex  eigen  wave  speed 


Suppose  that  a ^ W ^ b,  a and  b being  the  lower  and  upper 
bound  of  the  axial  velocity  respectively.  Equations  (3.28a, 
b)  imply 


0 < / (W  - a)(W 
'<1 


R 


- b)Qdr  = 


^ 2 

f W Qdr 


R. 


R. 


'2  "2 
(a  + b)  / WQdr  + ab/Qdr 


R 


1 


R 


1 


R 


{c  - (a  + b)c  + c^  + ab}  / Qdr 

^ R^ 


It  follows  that 


{c^  - l(a  + b)}^  + c?  < {l(a  - b)} 


(3.29) 
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This  is  equivalent  to  Howard's  semi-circle  theorem  for 
stratified  shear  flows,  saying  that  the  complex  wave 
speed  for  unstable  modes  must  lie  within  a certain  semi- 
circle, with  the  range  of  the  axial  velocity  for  diameter. 

It  was  Batchelor  and  Gill  (1962)  who  first  looked  at  the 
semi-circle  theorem  for  homogeneous  non-rotating  jet  flows. 
From  our  result  it  is  clear  theirs  can  be  extended  to  flows 
with  radial  density  variations. 

(b)  For  axisymmetric  perturbations  (m  = 0),  equations 
(3.27b.c)  reduce  to 

R2  R2 

/ WQdr  = Cj,  / Qdr  (3.30a) 


R 2 R2  ^2 

/W^Qdr  = (c^  + c?)  /Qdr  + / $|^l^rdr  (3.30b) 


If  a W £ b,  then  equations  (3.30)  imply 
0 < / (W  - a) («  - b)Qdr 


= (c^  - (a  + b)c,  + c?  + ab)  /Qdr  + fiii\^rdr 

r » 1 n D 


Consider  only  centri  fugal  ly  stable  profiles,  i.e.,  <I>  £ 0 . 


One  concludes  that 
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{c^  - |(a  + b)}^  + < {l(a  - b)}^  (3.31) 

From  examples  in  Sections  5.1.1(b)  and  5.2.3,  the  existence  of 
this  semi-circle  theorem,  first  derived  by  Leibovich  (1969) 
and  reducible  to  Howard  and  Gupta's  result  for  constant 
density,  suggests  that  even  for  axisymmetric  perturbations 
shear  instabilities  may  still  exist  when  the  Rayl ei gh-Synge 
condition  is  satisfied,  i.e.,  4'  ^ 0. 

(c)  For  azimuthally  periodic  perturbations  ( k = 0) 
and  D(p^fj)  = 0,  equations  (3.27b,c)  reduce  to 

R2  ^ R2^ 

f QQdr  = c / Qdr  (3.32a) 

R-,  '^R, 


(3.32b) 


where 

2 

Q = Pq[^|  D i|j|  ^ + i4^|  ^]r  > 0 
m 

c = ‘il 
m 

Let  a and  b be  the  upper  and  lower  bound  of  the  angular 
velocity  Q,  and  making  use  of  equations  (3.32)  we  obtain 
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For  centr i fugal ly  stable  profiles  with  Dp^  > 0 we  have 

{Cr  - ]-(a  + b)}^  + < {l(a  - b)}^  (3.33) 

Unfortunately  the  requirement  that  D(p^J^)  = 0 makes  this 
result  lack  generality. 

(d)  For  the  narrow  gap  approximation,  equations 
(3.27b,c)  take  the  form 


Rj 

/ (W  + s rf2)Qdr  = 


R 


1 


^2 

c^/Qdr 

^1 


(3.34a) 


/ (W  + s rj^)Qdr  = (c^  + c?)  / Qdr  + / G-^^^^-rdr  (3.34b) 


Ro  R 


2 I,,,  I 2 


^ R 


R^  k‘ 


Here  s = s 


m 


0 kR 


is  assumed  to  be  finite. 
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Letting 


a < W + < b 


then 


0 < 


/\(W  + s rf2)  - a][(W  + s^rfi)  - b]Qdr 
R2 


2 


(a  + b)c„  + + ab)  / Qdr  + /G — ^ 

R,  R. 


rdr 


For  centr i f ugal ly  stable  profiles  where  G > 0,  it  requires 


The  semi-circle  theorem  here  states  that  the  complex  eigen 
axial  wave  speed  c = ^ unstable  solutions  under  narrow 
gap  approximation  must  lie  inside  a semi-circle  which  has 
the  range  of  the  "swirling  velocity"  (W  + s^rQ)  as  diameter. 

Two  phenomena  can  be  observed  from  the  derivation 
of  the  semi-circle  theorem  for  the  four  special  profiles 
(and  for  two-dimensional  stratified  shear  flows  as  well) 
mentioned  above.  First,  it  appears  that  the  semi-circle 
theorem  yields  a bound  upon  the  complex  wave  velocity 
(either  c or  c)  only  for  shear  instabilities  with  centrifu- 
gally  stable  profiles;  second,  this  velocity-dependent 
bound  is  explicitly  independent  of  the  boundary  or  density 


- ^(a  + b)}^  + < {^(a  - b)}^ 


(3.35) 


effects . 
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3.3.2  The  Semi-circle  Bound  for  Rigid  Body  Rotation 

The  semi-circle  theorem  we  have  just  discussed  is 
the  bound  of  complex  wave  speeds  for  centri fugal ly  stable 
profiles  encountered  in  shear  motion.  A similar  criterion 
can  also  be  found  for  instability  of  rigid  body  rotation 
subject  to  azimuthally  periodic  disturbances. 

Considering  k = 0 and  in  equation  (2.11), 

one  obtains 


Dlp^rDc)))  - [— — + A(DPq)](()  = 0 


(3.36) 


with  the  boundary  conditions 


^(R^)  = cf)(R2)  = 0 


where 


(|)  = ru 


A 


2 


(3.37a) 


0 


with  the  imaginary  part 


A . 
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Equation  (3.36)  with  boundary  conditions  forms  a Sturm- 
Liouville  system  which  requires  that:  (1)  A be  real;  (2) 
A and  Dp^  be  of  opposite  sign.  The  first  condition  for 
c > 0 implies 


^0  \ 2 
2 ’ 


0x2 


(3.38) 


The  complex  azimuthal  wave  velocity  for  unstable  modes  must 
lie  exactly  on  a semi-circle  with  the  rotating  velocity  for 
diameter.  Meanwhile  (3.37a)  can  be  rearranged  as 


c ^ A - 1 ± /I  - A 

Q A 

0 


(3.39) 


c will  be  complex  whenever  A > 1,  so  that  Dp^  < 0 for  in- 
stability followed  from  the  second  condition  of  the  Sturm- 
Ljouville  system.  Note  that  if  Dp^  > 0 then  A <0  which  implies 
stability.  This  sufficiency  condition  for  stability  is  also 
obtainable  from  condition  (3.15)  by  considering  rigid  body 
rotation  for  azimuthal  periodic  disturbances. 


3 . 4 The  Growth  Rate 

The  semi-circle  bounds  we  have  discussed  are  the 
bounds  of  the  complex  wave  speeds  for  particular  types  of 
velocity  and  density  profiles  under  the  influence  of  shear 
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instability.  A general  upper  bound  on  the  growth  rates 
for  both  centrifugal  and  shear  instabilities  can  also  be 
derived  from  the  general  integral  equation  (3.12). 

Equation  (3.24c)  for  possible  instability  requires 

F|N|^  + 2G(kW  + mJl  - w^)  = 0 (3.40) 

at  least  at  one  point  in  (R-|,  R2)-  Equation  (3.40)  is  a 
quadratic  equation  for  the  complex  u)  and  can  be  rearranged 
into 

kW  + mfl  - 0)  = -G  ± </G^  - F^u)?  (3.41) 

r I 

The  left-hand  side  of  (3.41)  is  a real  quantity,  and  it 
follows  that  the  function  inside  the  radical  should  be 
positive  to  have  real  solutions,  i.e., 

r2  .2  2 ^ n 

G - F ^ 0 
or 

0) . < max  I f1  (3.42) 

provided  {F}  does  not  change  sign  in(R^,  R^)- 

We  can  observe  immediately  from  (3.42)  that  this 
bound  fails  to  predict  growth  rates  when  {F}  becomes 
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sufficiently  small  anywhere  inside  the  flow  domain.  For 
instance,  becomes  infinity  for  a potential  vortex  with 
a superposed  constant  axial  flow,  which  is  actually  stable 
against  all  disturbances  except  for  negative  density  gradi- 
ents. 

On  the  other  hand,  equation  (3.42)  can  reasonably 
predict  the  growth  rate  for  solid  body  rotation  subject 
to  azimuthally  periodic  perturbations.  The  upper  bound 
of  the  present  flow  predicted  in  this  case  is 
which  is  exactly  equivalent  to  the  result  obtained  from 
(3.38) . 

Another  second  upper  bound  on  the  growth  rate  is  ob- 
tainable from  the  general  integral  equation  by  setting 
C = ^.  For  instability,  equation  (3.14)  can  be  written 
as 

^2  1 1 2 ^2  * 

/ {p  k^lNl^  + [G  - Ip  E(DN)2]}-^rdr  = - f p E|D%|^rdr  s 0 

° ^ ° |Nr  Ri 

(3.43) 


The  inequality 

p^k^lNl^  + [G  - Ip^E(DN)^]  < 0 

must  hold  somewhere  inside  the  region  of  the  flow.  Since 
(jj?  < |N1^,  we  find 
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1 [DW  + s(4fi  + rDf2)l^ 
4 1 + s2 


P 


} 


(3.44) 


0 


This  is  reducible  to  Howard  and  Gupta's  result  for  constant 
density.  Several  observations  can  be  made  from  this  upper 
bound  immediately.  Axial  and  angular  velocity  gradients 
destabilize  the  flow  and  definitely  increase  the  growth 
rate,  while  the  centrifugal  effects  stabilize  the  flow 
and  reduce  the  growth  rate.  No  instabilities  can  occur 
whenever  the  centrifugal  stabilizing  effects  are  suffi- 
ciently strong  to  suppress  the  destabilizing  ones  of 
shearing  rates  of  both  the  axial  and  tangential  velocity 
components . 


CHAPTER  4 

INTERFACIAL  CONDITIONS 

From  the  general  criteria  discussed  in  Chapter  3, 
stable  swirling  flows  will  exist  whenever  X 1 4-  is  satis- 
fied. However,  violations  of  this  sufficiency  condition 
do  not  necessarily  imply  instability.  In  order  to  deter- 
mine if  the  flow  is  unstable  when  condition  (3.19b)  is 
violated,  it  is  necessary  to  solve  the  second  order  dif- 
ferential equation  (2.11)  or  the  two  corresponding  first 
order  differential  equation  (2.9)  and  (2.10)  exactly. 
Unfortunately,  examples  of  swirling  flows  for  which  the 
stability  problem  can  be  explicitly  solved  do  not  seem  to 
be  readily  found.  The  trouble  with  finding  analytical 
solutions  to  (2.11)  is  that  the  equation  possesses  a 
second  order  as  well  as  a first  order  singularity.  Equa- 
tion (2.11)  is  singular  whenever  N = kW  + mii  - w = 0 and 
only  series  solutions  can  be  obtained  by  the  Frobenius 
method  with  a restriction  that  the  swirling  velocity 
gradient  { DW  + srDfi} should  nowhere  vanish  inside  the  flow 
region.  A second  look  at  the  singularity  shows  the  dif- 
ficulty in  locating  the  possible  singular  points  which 
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are  in  terms  of  the  unknown  eigen  wave  speed. 

Since  explicit  solutions  to  the  eigenvalue  prob- 
lem, especially  for  unstable  modes,  are  difficult  to  find 
for  smoothly  varying  functions,  the  use  of  non-smoothly 
varying  profiles  can  be  sometimes  justified  as  an  approxi- 
mation to  related  smoothly  varying  ones.  The  analysis  of 
the  instability  for  broken-line  velocity  and  density  dis- 
tributions has  been  extensively  covered  in  the  literature 
and  can  be  employed  to  yield  many  instability  character- 
istics. In  general,  this  analysis  will  allow  simple  solu- 
tions to  the  governing  differential  equation  in  the  in- 
dividual region  which  possesses  continuous  si mpl e- vary i ng 
velocity  and  density  distributions.  We  will  proceed  to 
construct  the  stability  characteristics  of  such  broken- 
line  velocity  and  density  profiles. 

In  deriving  the  stability  equation,  we  have  not 
considered  the  possibility  of  the  existence  of  forces 
(such  as  surface  tension  and  centrifugal  effects)  arising 
from  the  possible  discontinuity  at  the  interfaces.  Modi- 
fications in  this  case  become  necessary  due  to  the  impor- 
tant roles  played  by  these  forces  which  do  not  exist  in 
continuous  velocity  and  density  distributions. 

Consider  surface  tension  T.  acting  at  the  inter- 

vJ 

face  j located  between  two  rings  of  fluids  at  r = Ry 
The  equation  of  motion  for  the  mean  flow  in  the  radial 
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direction  is 


dP  , 

0 nC 

-3 — - p rJ^ 
dr  ^0 


M 1 

5:  T.  {^)  6(r  - R.) 

j=l  J J ^ 


(4.1  ) 


where  6(r  - R.)  is  a Delta  function,  R.  is  the  position 
J J 

as  well  as  the  principal  radius  of  curvature  for  the  dis- 
turbed interface. 

Once  the  flow  is  disturbed,  the  principal  radii  of 
curvature  on  r - 0 and  r - z planes  are  specified  as 


2 


30 


(4.2) 


and  the  linearized  perturbation  equation  of  motion  in  the 
radial  direction  yields 


dP 


r + 1^-  i P„(d*u  - 2fiv)  - (p^  + p)rfi^ 


dr  3r  ^0 


X T,( 


1 1 


j = l j ^rQ  '^rz 


-)  6(r  - R.) 


(4.3) 


Subtracting  (4.1)  from  (4.3),  one  finds 
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+ pQ  - 2s^v)  - prn^ 


(4.4) 


In  the  absence  of  surface  tension,  equation  (4.4)  becomes 
identically  equal  to  equation  (2.5a).  The  validity  of 
the  equations  of  motion  in  the  9 and  z directions,  the 
continuity  equation  and  the  incompressibility  condition 
is  unaffected  by  the  existence  of  surface  tension  and 
needs  no  modification. 

Following  the  normal  mode  analysis  and  assuming 
infinitesimal  disturbances,  we  let 


at  the  interface  j.  Equation  (4.2)  can  then  be  linearized 
into  the  form 


r = R.  + n = R.  + n(R.j)e 

O J J 


i ( kz  + me  - o)t ) 


R 


re 


+ 


R 


rz 


2 2x  i(kz+me-u)t) 
m - K . ) n.e  ' 

u 0 


(4.5) 


where 


Assume  periodic  functions  for  all  the  perturbation  quanti- 
ties, i ^ e . , 

(0,  w,  P,  p)  = fu(r),  v(r),  w(r),  p ( r ) , p ( r ) ) e ’ ' 

(4.6) 

Equations  (4.4),  (4.5)  and  (2.5b  to  e)  can  be  transformed 
into  ordinary  differential  equations,  namely. 

Dp  + p [iNu  - 2av]  - prn^  = E -^  ( 1 - m - K^)nj6(r  - ) 


p + Po  ^[Nv  - iD*(rfl)u]  = 0 

p + ^[Nw  - i(DW)u]  = 0 (4.7) 

iNp  + (Dp^)u  = 0 

D*u+ik(sv+w)=0 

Two  first  order  differential  equations  can  be  deduced  from 
(4.7). 


(k^  + 2^)P 


TP, 


[k(DW)  + - D*]u 


r 


(4.8a) 
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NDp  + — p 


i { ( <I>  - PqN  ) u 


T . 

-i  n 

j = l 


M 

N E 


(^)  «(r  - R )) 

j 


(4.8b) 


Equations  (4.8)  are  exactly  equivalent  to  (2.9)  and  (2.10) 
except  for  the  surface  tension  term.  The  modifications  of 
the  stability  equation  due  to  surface  tension  arising  at 
interfaces  result  in 


N^[D(p^ED*)  - P^k^lu  + [NF  + G]u 


. M T . 

K^N  E (1 

j = i 


m^  - K^)(^)  6(r  - R.) 
J 


(4.9) 


Equation  (4.9)  without  surface  tension  reduces  to  (2.11), 
which  is  valid  at  points  not  on  a discontinuity  surface. 
Alterman's  (1961)  result  can  be  deduced  for  axi symmetri ca 1 
disturbances. 

Both  the  di spl acement  and  pressure  conditions  will 
be  discussed  mathematically  and  physically  as  follows: 


4.1  Kinematic  Interfacial  Condition 


Consider  two  rings  of  fluids  in  region  1 and  region 
2 with  an  interface  at  r = R.  Integrating  (4.8a)  from 
R - e to  R + t and  assuming  that  the  mean  value  exists 
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for  an  the  quantities  of  the  flow  and  that  their  de- 
rivatives should  be  bounded  at  the  interface,  we  obtain 

R+e  R+e 

lim  ^1  / [m(Dn)  + k(DW)]dr  - / (Du)dr  = 0 

e->0  ^ R R-e 

or 

^1  [m(^^2  - ^i)  + MW2  - W^)]  - [U2  - u^]  = 0 (4.10) 

R 

where 


N2^’ 


the  mean  value  at  the  interface.  Define 
<f>  = f(R,,)  - f(R_o) 


(4.11  ) 


a jump  which  a quantity  {f}  experiences  at  the  interface 
r = R.  Without  making  use  of  any  physical  arguments, 
equation  (4.10)  can  be  written  as 

<^>  = 0 '■7,)'’  (4.12) 

the  quantity  {^},  physically  interpreted  as  a perturbation 
displacement,  must  be  nowhere  discontinuous  at  the  interface 
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or  anywhere  inside  the  domain  of  the  flow. 

Condition  (4.12)  is  also  obtainable  by  the  follow- 
ing physical  argument: 

Let  the  radial  perturbation  displacement  be  n . 

The  analysis  of  the  normal  mode  method  allows  all  the  per- 
turbation quantities  to  admit  solutions  with  periodic  be- 
havior in  2,  0 and  t and  with  radially  dependent  ampli- 
tudes. Hence  we  assume 

n(r,  0,  z;  t)  = n(r)e'<^^  " ‘ (4.13) 

According  to  the  chain  rule 

dt  “ at  3r  dt  90  dt  9z  dt 

= K-O,  0 ma  + kW)n(r)e''"'4  " ‘ (4.14) 

= iNn(r)e’<^"  * 

By  definition,  the  total  derivative  of  the  radial  pertur- 
bation displacement  is  equal  to  the  radial  perturbation 
velocity  which  has  the  form 


u = u ( r ) e 


i(kz  + m0  - o)t) 


(4.15) 
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Equating  (4.14)  and  (4.15),  one  finds, 
u 

iN  (4.16) 

Equation  (4.12)  then  follows  from  the  physical  argument 
that  the  normal  component  of  the  perturbation  displacement 
must  be  continuous  across  the  material  fluid  interface  in- 
side the  flow  domain. 

4 . 2 Dynamic  Interfacial  Condition 

The  dynamic  boundary  condition  at  an  interface  be- 
tween two  fluids  has  been  carried  out  through  different 
arguments  for  continuous  and  discontinuous  velocity  and 
density  distributions.  The  application  of  an  integration 
technique  to  equations  governing  stability  across  inter- 
faces has  been  used  to  lead  to  explicit  interfacial  con- 
ditions for  problems  in  three-dimensional  parallel  flows 
(Chandrasekhar  1961)  and  two-dimensional  parallel  flows 
(Drazin  and  Howard  1966).  Alterman  (1961)  looked  at  the 
stability  characteristics  of  a heterogeneous  liquid  with 
radially  dependent  azimuthal  and  axial  velocities.  For 
antisymmetric  disturbances,  he  obtained  an  interfacial 
condition  for  two-region  flows  by  integrating  the  govern- 
ing second  order  differential  equation  across  the  inter- 
face between  the  two  fluids.  Kent  et  al . (1969)  used  the 
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same  method  to  obtain  a jump  condition  at  the  interface 
for  the  Kel vi n-Hel mho! tz  instability  around  an  isothermal, 
cylindrical  plasma  in  an  arbitrary  radially  dependent  den- 
sity and  electric  field.  Unfortunately,  this  application 
of  the  integration  technique  to  second  order  equations 
governing  stability  can  be  used  to  lead  to  simple  results 
for  only  some  particular  cases.  General  pressure  condi- 
tions cannot  easily  be  found  by  this  approach.  However, 
the  general  dynamic  interfacial  condition  can  in  fact  be 
derived,  using  the  integration  technique  across  the  inter- 
faces, from  one  of  the  two  first  order  differential  equa- 
tions leading  to  the  stability  equation.  Equation  (4.8b) 
can  be  rearranged  to  the  form 


Applying  the  integration  technique  across  the  interface 
at  r = R,  we  obtain 


Dp  + = HCDlP^rn^)  + - N^)](^) 

S -i  (1  - m^  - <5(r  - R.)} 

i = l R.  J " j 


(4.17) 


1 i m 

£->■0 


.(^)  6(r  - R)]}  dr  = 0 
R 
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or 

<p>  - i(^)  {<p  rn^>  - ( 1 - - K^)}  = 0 (4.18) 

IN  0 

< >,  defined  in  (4.11),  represents  a possible  jump  condi- 
tion at  the  interface. 

The  reason  why  the  general  dynamic  interfacial  con- 
dition is  difficult  to  obtain  by  integrating  the  second 
order  differential  equation  for  stability  has  now  been  re- 
vealed. Equations  (4.8a,b),  the  two  first  order  differ- 
ential equations  leading  to  the  stability  equation,  repre- 
sent the  kinemetic  and  dynamic  interfacial  conditions  re- 
spectively, which  couple  together  in  the  second  order 
stability  equation.  Except  in  special  cases  they  are  not 
likely,  in  general,  to  be  separated  if  integration  tech- 
niques are  considered  only  with  respect  to  the  equation 
governing  stability. 

Equation  (4.18)  may  also  be  reached  through  the 
following  physical  argument: 

Consider  two  regions  of  fluids  with  an  interface 
at  r = R inside  an  annular  ring  R-|  £.  r £ equation 

of  motion  (4.1)  is  valid  for  the  flow  in  both  regions. 

Let  the  subscripts  1 and  2 represent  the  quantities  in 
region  1 and  region  2 respectively.  Integrate  (4.1)  from 
the  inner  boundary  outwards  for  both  regions.  The  total 
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pressures  for  the  steady-state  are 


P-](r)  = / p-|  rJ2^dr  for  1 1 R 


R 


(4.19a) 


R-e  p R+e  j 

P„(r)  = lim  { / p,rf2fdr  - / -4,6  ( r - R)dr  + 
£->0  R,  ' ' R-e 


R + e 


2 

f p2rf^2^’^^ 


R 2 2 T 

/ p-|rfJ^dr  + / P2rn2clr  - for  R £.  1 R 


2 


(4.19b) 


where  T stands  for  the  surface  tension  evaluated  at  the 
interface.  From  (2.4)  and  (4.19),  the  total  pressures 
for  the  disturbed  state  evaluated  at  the  deformed  inter- 
face with  the  perturbation  displacement  n(r,  0,  z;  t)  can 
be  expressed  as 

R+n  2 ^ 

P-|(R  + n)  = / p-|rf2^dr  + p-j(R  + n) 


R 


(4.20a) 


+ Po ( R + n ) 


(4.20b) 
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In  the  perturbed  state,  the  total  pressures  evaluated  at 
the  deformed  interface  should  be  balanced  by  surface  ten- 
sion through  the  vortex  sheet  (see  Figure  4.1). 


Figure  4.1  Total  forces  acting  at  the 
deformed  interface. 


P,(R  + n)  = P„(R  + fi)  + T(^-  + 

' ^ '^re  '^rz 


(4.21) 


or 


R+n 

/ 

R 


p-|  rfi^  dr 


P^(R  +n) 


R+n 
/ P- 
R 


P2(R  + n) 


- T[i 


're 


rz 


-)] 


(4.22) 
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Since  both  velocity  and  density  profiles  are  bounded  and 
continuous  in  [R,  R + n]  and  the  displacement  n is  assumed 
to  be  small,  equation  (4.22)  can  be  written  as 

R + fj 

(P?  - Pi)  - T[k  - ^ + ^)]  = (Pnfi?  - PpS^o)  / rdr 
' R ^ ^re  '^rz  ' ' ^ ^ R R 

= 

= (4.23) 

Consider  (4.5),  (4.6)  and  (4.13)  for  periodic  solutions. 
Equation  (4.23)  can  then  be  simplified  to 

f 

<p>  + {<p^rfi^>  - -^(1  - m^  - K^)}  Tin  " 0 (4.24) 

0 r2  K 

Equation  (4.18)  follows  by  substituting  (4.16)  into  the 
foregoing  equation. 

The  physical  meanings  of  the  mathematical  steps 
from  (4.21)  to  (4.24)  are  in  fact  the  dissolution  of  the 
total  pressure  forces  acting  at  the  perturbed  interface 
into  the  force  components  acting  at  the  steady-state  inter- 
face, provided  the  perturbation  quantities  for  the  second 
or  higher  order  are  negligible.  This  procedure  can  be 
summarized  in  the  following  sketch: 
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PifR) 


( )n  — > 
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^ P2(R) 

\ 

/ 

-«r-/(P2Ri^2^n 


/ 


/ 


r = R 
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Figure  4.2  Dissolution  of  total  pressure  forces  into 
force  components  at  the  interface. 


In  Figure  4.2b,  we  consider  a differential  area  dA  = Rde 
at  the  interface  r = R.  The  centrifugal  forces  due  to 
the  rotation  of  the  corresponding  differential  fluid  vol- 
ume  ri(Rd0)  with  respect  to  the  axis  are  { RJ^^  ) (nRdO ) and 
{p2Rfl2)  fo*"  regions  1 and  2 respectively.  Summing 

up  all  the  radial  forces  acting  at  the  fluid  element,  one 
obtains  the  equilibrium  condition 


(P.|  + p^)(Rde)  + p^RJ^^(nRd0)  = {P2  + P2)(Rd0) 

+ p«Rfi?  ( fjRdO  ) + T(p-^ — + p-^  )(Rd0) 
2 2 


(4.25) 
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In  the  steady-state 

P^(Rde)  = P^CRde)  + ^(Rde)  (4.26) 

Subtracting  (4.26)  from  (4.25),  we  find 

(Pp  - Pi)  + [(PpRfip  - P-|Rf^?)n]  + T[^  ■ ^ R^^^  " ° 

^ ' R ^ ^ ' ' R re  rz 

This  is  exactly  equation  (4.23),  and  equation  (4.24)  follows 
for  periodic  solutions. 

It  should  be  emphasized  that  both  kinematic  and 
dynamic  interfacial  conditions  were  derived  from  the  same 
set  of  equations  for  heterogeneous  swirling  flows  without 
any  other  assumptions,  and  they  must  be  satisfied  when 
matching  techniques  for  broken-line  profiles  are  applied  at 
i nterf aces . 

Michalke  and  Timme  (1967)  examined  the  instability 
of  a potential  vortex  enclosing  a stagnant  fluid  core  for 
a homogeneous  density  system.  In  matching  the  interfacial 
condition  at  the  cylindrical  vortex  sheet,  they  only  con- 
sidered equating  the  perturbation  pressures  in  spite  of 
the  angular  velocity  jump  governed  by  equation  (4.18). 

Uberoi  et  al.  (1972)  pointed  out  the  error  by  taking  into 


account  the  first  order  perturbation  of  equilibrium  pressure 
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across  the  deformed  surface.  Their  argument  is  essentially 
the  following: 

Equations  (2.1)  and  (2.4)  may  be  written  as 


P . 
J 


/ p . rJ2  .dr  + p . 


1 


(4.27a) 


or 


j = 1,  2 


(4.27b) 


We  use  the  Taylor  expansion  with  respect  to  the  perturba- 
tion displacement  fi  and  neglect  the  second  or  higher  order 
perturbation  quantities.  Equations  (4.27)  then  yield 

dP.  r p ^ o ^ 

Pj  + ^Tr  ""  ^ p^rs^^dr  + Pj  + (Pj>"f^j)n  j - 1,  2 (4.28) 

■^1 

The  total  pressure  together  with  the  surface  tension  should 
be  balanced  at  r = R. 

R 2 ^ 2 ^ P ? 

/ pirJ^^dr  + p-|(R)  + [(p-|rf^^)n]  = / p^rQ^dr  + P2(R)  + 

R 1 R R 1 

+ T(^  + (4.29) 

R r0  rz 
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In  the  steady-state 
R R 

/ + T(^)  (4.30) 

Rl 

Subtracting  (4.30)  from  (4,29)  and  assuming  period  solu- 
tions, we  obtain  the  same  dynamic  interfacial  condition 
(4.18). 

For  piecewise  continuous  profiles,  one  can  join 
up  the  solution  in  the  proper  region  by  use  of  conditions 
(4.12)  and  (4.18)  at  interfaces.  This  method,  first  used 
successfully  by  Kelvin,  has  led  to  explicit  eigenvalue 
relations  for  many  problems,  which  serve  as  an  approxi- 
mation to  smoothly  varying  profiles  for  small  wave  numbers. 


CHAPTER  5 


SOME  ANALYTICAL  SOLUTIONS 

As  mentioned  in  Chapter  4,  the  difficulty  in 
obtaining  exact  solutions  to  equation  (2.11)  arises  from 
the  singularity  existing  in  the  equation.  The  location 
of  possible  singularities  and  solutions  in  terms  of 
known  functions  do  not  seem  to  exist.  Series  solutions 
by  the  Frobenius  method  are  only  obtainable  when  {W  + sr^  - w} 
and  the  swirling  velocity  gradients  {DW  + srDfi}  do  not 
vanish  simultaneously  anywhere  within  the  flow  domain. 

This  difficulty  can  only  be  overcome  either  by  getting  rid 
of  the  singularity  in  equation  (2.11)  (i.e.,  by  considering  a 
vortex  with  constant  density  and  zero  axial  flow  gradient) 
or  by  setting  N equal  to  a constant,  which  requires  constant 
swirling  velocity  profiles.  All  analytical  solutions  ever 
found  have  fallen  into  one  or  the  other  of  these  two  cate- 
gories: the  former  was  first  examined  by  Ponstein  (1959) 

for  a potential  vortex  superposing  a liquid  jet  surrounded 
by  air,  and  subsequently  by  Michalke  and  Timme  (1967)  for 
a vortex  enclosing  a stagnant  core  for  constant  density; 
the  latter  was  studied  by  Alterman  (1961)  for  the  capillary 
instability  of  a liquid  jet  streaming  within  a surrounding 
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fluid,  and  by  Reshotko  and  Monnin  (1965)  for  a two-fluid 
wheel  flow  in  the  absence  of  surface  tension.  Solutions 
not  falling  into  these  two  categories  have  not  been  found. 

We  now  present  examples  of  exact  solutions  of 
certain  swirling  flows  with  discontinuous  velocity  and 
density  profiles  subject  to  arbitrary  infinitesimal  per- 
turbati ons . 

5 . 1 Solutions  in  Terms  of  Step  Approximations  to  Continuous 
Profiles 


Consider  one  ring  of  fluid  with  constant  density 

p.,  constant  angular  velocity  Q.  and  constant  axial  veloc- 
J J 

ity  Wj  in  region  j.  Equations  (2.9)  and  (2.10)  for  the 
profiles  under  consideration  can  be  simplified  to 


- njD*]u  = (k^  + ~)p  (5.1) 

- n^)u  = (n^Dp  + ^p)  (5.2) 


Eliminating  u from  the  above  two  equations,  one  obtains  a 
differential  equation  for  the  fluctuation  pressure 

2 

D^p  + iop  - [k^q^  + !^]p  - 0 (5.3) 

r 


where 


~T 


n . 
J 


J 


+ m - 


U) 

“j 


1 
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Equation  (5.3)  is  the  modified  Bessel  equation  with 
general  solutions  in  terms  of  the  linear  combination 
of  the  modified  Bessel  functions  of  the  first  and 
second  kinds,  i . e . , 

P = ^ ^ ^ ^ (5.4) 


here 


>=j  ' "j'''' 

The  normalizing  factors  k"’’^  and  k"'  are  here  introduced 
in  order  to  allow  k ^ 0 for  m ^ 0,  which  gives 


lim  p = -ip.S^^(n^  - 4){A  r^  + B r"*^} 
k->0  J J J 


where 


* 

A 


m 


= A- 


0^^  I 

Z m ! 


Substituting  (5.4)  into  (5.2)  leads  to  a solution  for  per- 
turbation velocity  of  the  form 


-m 

u = kq.N.iAk  "’[-J-J2 — J 

J s) 


. + 


/X 


m Xi^  (xJ 

Rk^r  J ^ J . + 
L K (x.) 
m'  '^j  ' 


"j 


(5.5) 
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To  approximate  continuous  velocity  and  density 
profiles,  consider  a multiple  layer  system  with  M rings 
of  fluid  and  M-1  interfaces  inside  the  annulus  between 
two  cylindrical  walls  of  radius  and  R|y|,  the  velocity 
and  density  distributions  being  step  functions  as  shown 
in  the  following  sketch: 

Pq,  w,  n 


At  the  interface  r = R.,  two  rings  of  fluid  ,j  and  j + 1 
are  connected.  The  displacement  condition  (4.12)  with 
solutions  (5.5)  in  region  j and  j + 1 gives 
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mjj  J m'^jj'j 


+ — ]K  (q.^.K.)  = 0 
nj+i''  m'^j+1  y 


(5.6) 


while  the  pressure  condition  (including  surface  tension 
T.  acting  at  the  interface)  yields 

\J 


j 


«jk">{(„2  - 4)  . , r (1  - „2  - <^)]iyq  K ) 

J m j'  ^ J J J 


m'^j+1  j' 


■a . 


(5.7) 
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where 


r . 

j 


T . 


2 3 


a 


j + 1 , j 


2 


Equations  (5.6)  and  (5.7)  represent  2M  linear 
homogeneous  equations  in  the  variables  of  interfacial 
displacements  and  pressures.  Be  reminded  that  the  fluctu- 
ation velocity  should  vanish  at  the  walls;  hence 


A,k-"'{^  + 
1 n, 


q,  K I ' (q,  K ) 
^1  0 m ^1  0 


m 


(q^  ^ ^m'^ro'  ' "1  n 


) + + 


q,K  K (qiK  ) 
^1  0 m'^1  0 


o' 


}K  (q,K  ) - 0 

m ^1  o' 


(5.8) 


A„k-"'{Sl 
M "m 


KWtAi  M n. 


'M 


q.,KwK ' (q^K..,) 
M m M 


}K  (qwKM)  = 0 
m'^M  M' 


(5.9) 

For  nontrivial  solutions,  the  determinant  of  the  coeffi- 
cients of  (5.6)  to  (5.9)  must  be  zero.  This  yields  the 
di spers ion  relation 
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1,0 

“^l  ,0 

1,1 

'^1,1  ‘^2,1 

■'^2,1 

* 

* * 

★ 

1,1 

*^1,1  '^2,1 

■'^2,1 

I.  . K.  . -I.^T  . -K.^t  . 

J,J  J,J  J+1  .J  J+1  .J 

★ ★ ★ ★ 

I.  . K.  . -I.^T  . -K.^t  , 

J,J  J,J  J+1,J  J+1 


(5.10) 

where 


q.K.K' (q -K.)  p 


I = r?ltf , X 
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I*  . = {(n?  - 4)  + + r.(l  - - K^)]}I  (q.K.) 


J,J 


•n.  TJVV""  J 

J rn  J J 


m'  'J  J' 


K*  . = ((n^  - 4).  - qC  - 


J,J  J 


I*  =a  f(n2  -4)  + r?ni  + 

‘j+l.j  “J+I,j‘'"j+1  ' 'nj^, 


J*  ,.j i -ll  1 (n  y'  1 


4)  + [Jm  + ( ) 


(5.11) 

The  generalized  equations  which  describe  determinant  (5.10) 
are 


-m 


m 


m. 


m„ 


¥ ¥j^¥¥j-¥i^¥^ij-¥i'V 


(5.13) 


j = 1,2,  M-1 


with  boundary  conditions 


^1*^  ^1 ,0  ^1*^  “^l  ,0  ° 


Awk'^^lM  M + m = 0 

M M,M  M M,M 


(5.14) 
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Determinant  (5.10)  is  of  order  2(M  + 1)  with  the 
non-zero  elements  appearing  in  a systematic  form  around 
the  diagonal.  Equations  (5.12),  (5.13)  and  (5.14)  give 
the  general  solution  for  discrete  strata  for  any  number 
of  rings  of  fluid.  Continuous  stratifications  can  be  ap- 
proximated by  letting  the  number  of  rings  approach  infinity 
(M->“),  with  the  thickness  of  each  ring  of  fluid  approaching 
zero.  However,  one  should  keep  in  mind  that  the  velocity 
jump  at  interfaces  for  step-wise  approximations  of  continu- 
ous functions  will  yield  the  Kelvin-Helmholtz  instabilities 
for  large  wave  numbers,  while  the  corresponding  continuous 
distributions  are  quite  stable  for  all  modes.  This  dif- 
ference arises  from  the  strong  shearing  rates  existing  in 
the  sharp  shear  layer  created  at  the  vortex  sheet.  Instabil- 
ities produced  by  the  Kelvin-Helmholtz  effects  due  to  the 
relative  motion  at  interfaces  are  too  strong  to  be  stabil- 
ized for  short  wavelengths.  Condition  (3.19)  would  be  vio- 
lated at  certain  points  for  the  discrete  strata  while  it  is 
actually  satisfied  for  the  corresponding  continuous  profiles. 
This  phenomenon  is  of  significant  importance  especially 
when  the  wave  numbers  are  sufficiently  large.  The  approxi- 
mation of  continuous  distributions  by  many  rings  of  fluid 
offers  a convenient  approach  to  eigen  solutions  for  arbi- 
trary continuous  profiles  of  which  exact  solutions  are  not 
easily  obtained.  With  the  use  of  electronic  computers. 
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reasonably  good  approximations  for  eigenvalue  w = + ico^ 

are  expected  when  the  thickness  of  each  layer  is  much  less 
than  wavelengths  of  the  distributions  under  consideration. 

The  order  of  secular  determinant  (5.10)  is  re- 
ducible somewhat  by  a useful  theorem  (Pipes  1 970).  Con- 
sidering K.J  Q and  as  the  two  pivotal  elements  for 

the  first  and  second  reductions  of  the  determinant,  one 
obtai ns 
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1,1 


1.1 


J,J  J+1.J  J+l.J 


* * * 
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■j+1  ,j 


M-1,M-1  “^M-^M-l 


★ ★ 


M-1  ,M-1  "^M-l  ,M-1 


Sl.M-l  " ^M^M,M-1 


(5.15) 
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M ,f1 


(5.16) 
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If  the  inner  cylinder  vanishes  q will  diverge  to 
infinity  while  I-|  q remains  finite,  with  the  result  that 
will  approach  zero.  For  r^ , will  vanish  due  to 

the  divergence  of  at  infinity.  The  order  of  the 

determinant  has  now  been  reduced  to  2M,  where  M repre- 
sents the  number  of  rings  of  fluid  subdivided  between 
two  solid  cylinders.  For  some  finite  numbers  of  M,  one 
can  obtain  the  secular  relation  by  following  the  same 
track  to  keep  on  reducing  the  order  of  the  determinant. 

With  M = 2 for  two  region  flows  with  an  inter- 
face at  r = R,  the  secular  relation  reduces  to  a two  by 
two  determinant 


1 ,1 

'^2,1 

^2^2,1 

•k 

★ 

★ 

★ 

o 

II 

1 ,1 

- 

'^2,1 

^2^2,1 

(5.17) 


Expand  the  determinant  (5.17)  by  (5.11).  The  secular 
relation,  after  some  algebraic  procedures,  can  be  written 


as 


(n 


^ ) 


4)[I  (q,K) 
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where 


= a3^  - 1 - r,  (1  - - K^) 


(5.18) 
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q,  K,  I ' { q,  K,  ) 
2m  , 1 m' ^1  1 ' 


HI  I I n T / 
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r2m  , ^2^2’^m^^2^2^.,  , . 

L—  + v-~r-^-  j TJK^(q2'^2^ 


^ r2m  . 
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a = 


K = kR 


2 3 


'Ij  ' 7'  - \ "j  = - f: 

rij  J j 


W . 


j = 1 ,2 


Here  we  consider  R.|  as  the  inner  boundary  instead  of  R^, 
and  R represents  the  position  of  the  interface. 

Since  the  complex  eigen  amplification  factor  falls 
inside  the  arguments  of  the  Bessel  functions,  general  fea- 
tures of  the  system  are  difficult  to  see  without  the  use 
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of  a computer.  However,  some  important  results  can  be 
obtained  by  using  an  analytical  expansion  for  suffi- 
ciently large  or  small  axial  wave  numbers.  In  those 
cases,  the  Bessel  functions  reduce  to  an  exponential 
form  or  a power  series.  Stability  boundaries  and  growth 
rates  for  several  special  cases  are  presented  in  the  fol- 
lowing sections. 

5.1.1  Axisymmetric  Disturbances  (m  = 0) 

For  m = 0 equation  (5.18)  is  simplified  into 


(n;  - 4)[1  - 


("2  - 4)[1 
^2'-K„(q2K) 


K^(q2'<2)  io(q2-^) 

Ko(q2^2^  i;(q2<). 


= ic[a6^  - 1 - T-^{^  - K^)]  (5.19) 


Since 


r(z) 


Il(z) 


k;(z)  =-k^(z) 


Equation  (5.19)  can  further  be  reduced  to  the  form 


q^F^N^  + aq2F2N2  = - 1 "fid  - k^)] 


(5.20) 
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wi  th 


Il(qi<l)  Kp(q^K) 

( q] *^1 ) ^ 0 ^ 1 

1 1 1 ( q 1 ^ 1 ( q 1 ) K-|  { q 1 k) 

1^1  ( q 1 1^1  ) 1 1 ( q 1 


^0  ^ *^2  ^ 1 ( q 2 1^2  ^ 

2 K-j  ( q 2 ^ ) ^1  ( q 2 ) K-j  (CI2K2) 

1^1  ( q 2 ^ ^ 1 ( q 2 *^2  ^ 


The  influence  of  axial  velocities  on  the  stability  charac- 
teristics for  a xi symmetr i c perturbations  will  be  discussed 
as  follows: 

r 

( a ) Constant  axial  velocities 

. 'T  ' 

If  we  look  for  solutions  with  no  axial  flow  gradi- 
ent, i.e.,  = W^,  equation  (5.20)  takes  a simple  form 


^0 


w,) 


2 RJl^[a8^ 


- 1 - 


Fid  - 


K^)] 


c = 


cO 

k 


k ( q 1 F 1 ctq  2 ^ 2 ^ 


(5.21) 


Since  the  real  parts  of  the  modified  Bessel  functions  of 
the  first  and  second  kinds  are  monotonical ly  increasing 
or  decreasing  functions  respectively,  it  follows  that  the 
real  parts  of  both  and  are  positive  definite  for 
£ K £ <2*  Therefore,  the  flow  is  stable  whenever  the 
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term  inside  the  square  bracket  is  positive,  i.e., 

a6^  - 1 - r^(l  - ic^)  > 0 (5.22) 

This  result  was  first  obtained  by  Alterman  (1961)  with 
neither  the  inner  nor  the  outer  boundary.  It  is  not 
surprising  that  the  solid  walls  have  no  effects  on  the 
stability  boundary.  In  fact,  the  flow  we  are  now  con- 
sidering falls  into  the  category  of  centrifugal  instabil- 
ity, in  which  the  stability  characteristics  do  not  depend 
on  either  the  walls  of  the  cylinders  or  the  axial  wave 
numbers,  when  surface  tension  is  neglected.  This  is  analo- 
gous to  the  Rayl ei gh-Tay 1 or  instability  in  a gravitational 
force  field;  a statically  unstable  profile  is  going  to 
initiate  instabilities  despite  the  presence  of  boundaries 
or  the  effects  of  varying  wave  numbers.  If  we  look  back 
to  the  derivation  of  the  Rayleigh-Synge  criterion  by  dis- 
cussing the  characteristics  of  the  partial  differential 
equations  in  the  beginning  of  Chapter  3,  we  realize  that 
the  solid  boundaries  nowhere  affected  the  criterion  ex- 
plicitly and  no  assumptions  were  made  on  the  form  of  the 
perturbation  stream  function  either  (i.e.,  independent 
of  the  axial  wave  numbers). 

Instead  of  solving  the  governing  stability  equa- 
tion, the  stability  boundary  for  axisymmetric  disturbances 
in  the  absence  of  surface  tension  and  axial  flow  gradients 
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can  be  alternatively  located  through  a simple  integra- 
tion technique.  Since  both  flows  are  stable  in  their 
individual  regions  1 and  2,  the  instabilities  governed 
by  the  Rayl eigh-Synge  criterion  can  only  be  initiated 
right  at  the  interface.  In  applying  the  integration 
technique  with  respect  to  this  criterion  across  the 
interface,  one  finds  the  stability  condition  is 

lim  / D[p  (r^n)^]  = > 0 

e-^0  R-e  ° ^ ^ ' ' 

i . e . , 

aB^  - 1 > 0 (5.23) 

Both  (5.22)  and  (5.23)  are  exactly  equivalent  except  for 
the  surface  tension  term.  This  method  can  be  used  to 
locate  stability  boundaries  for  constant  axial  velocity 
subject  to  axisymmetric  perturbations  for  any  two  region 
flows  provided  the  flow  in  each  individual  region  is 
stable.  Stability  boundaries  are  shown  in  Figure  5.1. 
Surface  tension  stabilizes  or  destabilizes  the  flow  de- 
pending on  whether  ic  is  greater  or  less  than  one. 

( b ) Effects  of  axial  velocity  gradients 
As  previously  mentioned  in  Chapter  3,  the  shear 
effect  of  axial  velocity  can  cause  instabilities  of  the 
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n 


without  axial  velocity  gradients. 
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flow  even  when  the  Rayl eigh-Synge  condition  is  satisfied. 
A special  case  from  (5.19)  will  serve  as  a demonstration 
to  this  phenomenon . 

For  an  inner  boundary  at  r = 0 and  outer  boundary 
at  r = “ , equation  (5.19)  without  surface  tension  takes 
the  form 


2 2 
n.|  - 4 2 *^2  " ^ 

q-1  iglq-i <1  ~ q^i^^Tq^ 

“ITT^ 


(5.24) 


Consider  the  neutrally  stable  profile  (in  the  centrifugal 
2 

sense)  a6  = 1 and  expand  the  modified  Bessel  functions 
in  their  asymptotic  forms  valid  for  sufficiently  large  k. 
Equation  (5.24)  then  yields 


4 + 


4 


0 


(5.25) 


Squaring  (5.25)  and  rearranging  it,  one  finds 

[(W^  - c)2  - a(W2  - c)2][(w^  - c)2  + a(W2  - cf-  4(^)^]  ^ 0 


For  unstable  solutions,  it  requires 

[(W^  - c)2  + a(W2  - c)^  - 4(^)2]  = 0 


(5.26) 
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or 


c 


k ( 1 + a) 


(5.27a) 


where 


k(W2  - W^) 


Unstable  modes  are  expected  whenever 


(5.27b) 


This  example  clearly  demonstrates  that  the  Rayleigh- 


Synge  criterion  fails  to  be  valid  if  the  axial  velocity 
gradient  is  taken  into  consideration  (Chandrasekhar  1 960b), 
and  that  the  shear  instability  for  centri f ugal ly  stable 
profiles  is  governed  by  condition  (3.19). 


growth  rate  and  the  wave  speed.  It  is  interesting  to  note 
that  the  complex  wave  velocity  falls  inside  a semi-circle 
as  predicted  in  category  (b).  Section  (3.3.1)  of  Chapter 
3.  Furthermore,  all  unstable  solutions  for  short  waves 
lie  exactly  on  a smaller  semi-circle  governed  by  the  equa- 
ti  on 


Figures  5.2a,b  show  the  stability  boundary,  the 


I 2 

i 'max 


(5.28) 
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Density  effect  on  the  complex  wave  speed  for 
the  centrifugally  stable  profile  = 1 
= 1). 


Figure  5.2b 
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where 

I I I I 

^ ^ i Lax  2k  L-|  ^ 

obtained  by  differentiating  the  imaginary  part  of  (5.27a) 
with  respect  to  a for  unstable  solutions.  Here  the  plus 
or  minus  sign  in  (5.28)  depends  on  whether  the  value 
(W.|  - is  positive  or  negative  respectively.  Either 

large  or  small  density  ratio  has  a stabilizing  effect. 
Fora<  0.19  or  a = °°  , all  the  shear  instabilities  will 
be  completely  stabilized  by  this  density  effect  and  the 
flow  will  be  stable.  A liquid  jet  directed  into  a rotat- 
ing air  environment  or  a rotating  air  jet  blowing  into  a 
heavy  density  medium  is  stable  under  axisymmetric  dis- 
turbances provided  surface  tension  is  neglected.  Growth 
rates  are  expected  to  increase  as  the  densities  of  two 
fluids  approach  each  other,  and  to  reach  a maximum  at 
a = (y^  + 4)/(y^  - 4)  = 1 .38. 

5.1.2  Azimuthally  periodic  perturbations  (k  ^ 0) 

From  the  general  discussions  of  the  criteria  for 
rotating  flows  in  Chapter  3,  it  is  suggested  that  the 
two-dimensional  perturbations  for  k = 0 grow  more  rapidly 
than  the  corresponding  three-dimensional  ones.  Azi- 
muthally periodic  disturbances  in  rotating  flows  can 
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therefore  not  be  neglected  in  stability  investigations. 

For  small  axial  wave  numbers  (k  <<  1),  the  modi- 
fied Bessel  function  in  (5.18)  can  be  approximated  by 
their  algebraic  forms  for  small  arguments.  The  secular 
relation  for  the  two-region  flow  is  then  rearranged  into 
the  following  form  after  some  algebraic  procedures. 

n.|(f.|n.|  - 2)  + a3^n2(f2n2  2)  = m[a3^  - 1 - r.|(l  - m^  - k^)]  (5.29) 


where 


1 + 


1 - 


> 1 


e 


1 


j - 1,2 


Equation  (5.29)  is  a quadratic  form  for  w with  the  solution 


CO 

^2, 


W 

1]  + ot[f2(k—  + mB)  + B] 


± /a 


f,  + afj 


(5.30) 


where 

A = (f.|  + af2)[(:p — m)  + aB^(:p  + m)  - mr^(l  - m^  - k^)] 
- af^f^Ey  + (:^  - m)  + B(^  + m)]^ 
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The  stability  boundary  is  described  explicitly  by  the 
characteristic  equation 

A = 0 (5.31  ) 

For  unbounded  flow  (R-j->0  and  R2^“)  both  boundary 
factors  f.|  = f 2 = 1 » and  (5.30)  reduces  to 


0) 


W,  w. 

(kj^  + m - 1)  + ot(k^  + m3  +3)  ± /A 

1 + a 


(5.32a) 


with  the  stability  boundary  at 


A = (1  + a)[(l  - m)  + a3^(l  + m)  - mT.j(l  - - k^)] 


- a[y  + (1  - m)  + 3(1  + m)]  = 0 


(5.32b) 


We  make  several  observations  from  (5.32b)  and  summarize 
them  as  follows: 

( a ) Effects  of  azimuthal  wave  numbers 
For  k = 0 and  in  the  absence  of  surface  tension, 
equation  (5.32b)  can  be  solved  for 

cx  ( m ^ - 1 ) ±Jma{m^  - 1 ) (a^  - ^) 
a (iT  + 1 ) ( m - a) 


3 


(5.33) 
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Stability  boundaries  for  different  values  of  m are 
plotted  out  and  compared  to  that  of  axisymmetric  dis- 
turbances in  (5.23).  These  axisymmetric  disturbances 

2 

will  be  stable  whenever  a6  >;  1 and  unstable  when 
2 

aB  <1.  At  the  same  time,  the  stability  boundaries 
corresponding  to  azimuthal  modes  are  observed  to  lie 
partially  in  the  region  that  is  stable  to  axisymmetric 
modes.  This  means  that  the  present  flow  can  be  quite 
stable  to  axisymmetric  disturbances  in  part  of  the  a-B 
plane  and  yet  unstable  to  non-ax i symmetr i c modes.  Such 
a centr i f ugal ly  stable  profile  can  be  upset  by  the  shear 
effects,  once  they  are  sufficiently  strong.  Indeed,  as 
m->-“  such  azimuthally  periodic  modes  will  amplify  in  time 
everywhere  except  along  the  line  B = l5  a^l.  These 
results  clearly  demonstrate  that  the  Rayl ei gh-Synge 
criterion  is  generally  neither  a necessary  nor  a suffi- 
cient condition  for  the  stability  of  a heterogeneous 
swirling  flow  when  allowance  is  made  for  non-ax i sym- 
metric disturbances,  and  show  that  the  non-axis ymmetric 
modes  of  instability,  especially  corresponding  to  k = 0 
which  amplify  most  strongly  in  this  case,  cannot  be  ne- 
glected when  discussing  the  stability  of  such  flows. 

Figure  5.4  shows  the  instability  patterns  for 
different  azimuthal  periodic  modes  when  B = 1.  The 
stream  function,  which  represent  the  paths  of  moving 
particles,  has  the  form 


S^(r) 
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Figure  5.3  Stability  boundaries  for  axisymmetric  (m=0)  and  azimuthal 1y  periodic  (k=0) 
perturbations . 
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Figure  5.4  Stability  patterns  of  rigid  body  rotation  for  azimuthal ly  periodic  modes 
with  a density  difference  at  the  interface. 
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for  r < R 

for  r > R 

These  patterns  would  be  observed  in  a reference  frame 
rotating  at  the  constant  angular  velocity  of  the  flow. 

( b ) Effects  of  surface  tension 

Equation  {5.32b)  for  k = 0 can  be  rearranged  as 

a(m^  - 1 ) ± fjam{a  + l)(m^  - l)[(a-l)  - (a  - m)(m  + 1)1’^] 

^ " a(m  + 1 ) (m  - a) 

(5.34) 

Figure  5.5  shows  the  effects  of  surface  tension  on  the 
stability  boundaries  for  different  values  of  m in  the 
a-  B plane.  Several  observations  follow  immediately 
from  the  behavior  of  the  curves.  Surface  tension  sta- 
bilizes the  flow  especially  when  the  relative  velocity 
at  the  interface  is  small.  This  stabilizing  effect  will 
not  be  significant  when  stabilities  are  dominated  by 
strong  shearing  rates  existing  at  the  interface,  i.e., 
two  rings  of  fluids  rotating  in  opposite  direction  or  one 
ring  rotating  a lot  faster  than  another.  The  effect  of 
surface  tension  becomes  significantly  important  as  rn 
increases.  The  stable  region  in  Figure  5.5  will  expand 
more  and  more  as  m becomes  larger  and  larger  and  eventually 


■ (^)*^  cosmS 
R ^ 

(^)  cosme 
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Figure  5.5  Effects  of  surface  tension  on  stability  boundaries  for 
azimuthally  periodic  disturbances. 
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it  will  cover  the  whole  a- 6 plane  as  . This  means 

that  the  most  unstable  mode  cannot  exist  under  the 
influence  of  surface  tension.  The  capillary  forces 
existing  at  the  interface  are  sufficiently  strong  to 
damp  down  the  onset  of  instability  corresponding  to 
very  short  azimuthal  wave  lengths. 


( c ) Growth  rates  corresponding  to  the  complex 
semi -ci rcl e 

Consider  a > 1 and  a3  = 1 for  azimuthally  periodic 
disturbances.  Equation(5.32a)  without  surface  tension  can 
be  simplified  into 


with 


1 + a 


■^Vm  ( g - 1 ) ^ - 1) 

/mg( 1 + g) 


(5.35a) 


The  flow  will  be  unstable  whenever 


g 


m + 1 
m - 1 


i . e . , 


e < 


m - 1 
m + 1 


(5.35b) 


The  solutions  for  complex  wave  speed  are  shown 
in  Figure  5.6.  As  previously  predicted  by  (3.33),  the 
complex  azimuthal  wave  velocity  {t  = ^)  lies  within  a 
certain  semi-circle  of  diameter  equal  to  the  range  of 
rotating  velocity.  The  growth  rate  ampl i fies as  m increases 


101 


Figure  5.6  Effect  of  azimuthal  wave  numbers  on  the  complex 
eigen  phase  velocity  for  azimuthal ly  periodic 
disturbances. 
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until  it  touches  the  periphery  of  the  semi-circle  for 
m-*-",  which  represents  the  upper  bound  upon  the  maximum 
growth  rate  for  the  present  flow. 

(d ) Effects  of  axial  flows  for  small  k 

The  presence  of  axial  flows  is  irrelevant  for 
k = 0.  These  effects  are  worth  considering  however  for 
small  axial  wave  numbers  (for  example,  k = 0.1)  if  vari- 
ations of  axial  flow  are  significant.  Equation  (5.32b) 
without  surface  tension  yields 

(1  + a)[(l  - m)  + ag^(l  + m)]  - a[y  + (1  - m)  + g(l  + m)]^  = 0 

(5.36) 

From  Figure  5.7  we  will  notice  immediately  that  the  stabil- 
ity boundaries  depart  from  the  point  a = 1,  3 = 1 once  the 
axial  flow  difference  is  taken  into  consideration.  The 
lower  the  azimuthal  wave  number,  the  more  the  stability 
boundaries  depart  from  that  point.  Another  phenomenon  is 
that  the  positive  y-|  has  a destabilizing  effect  when  two 
rings  of  fluids  are  rotating  in  the  same  direction  and 
stabilizing  effects  when  they  are  rotating  in  opposite 
directions.  These  effects  will  diminish  for  larger  azi- 
muthal wave  numbers  and  finally  have  no  influence  on  the 
stability  boundaries  as  m 
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small  axial  wave  numbers. 
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( e ) Effects  of  boundaries 
It  has  been  shown  previously  that  the  solid 
walls  do  not  have  any  effect  on  the  stability  boundar- 
ies if  only  axi symmetri c disturbances  are  considered, 
since  the  centrifugal  instability  depends  only  on  the 
centrifugal  pressure  balance  between  two  adjacent  fluid 
elements.  However,  this  argument  will  not  hold  any 
longer  once  shear  effects  become  a dominant  factor  in 
the  initiation  of  instabilities. 

The  stability  boundary  reduced  from  (5.31)  in 
the  absence  of  surface  tension  becomes 


a(mfi  - l)(mf2  + 1)  ± /am(ot  - l)(f^  + af2)(mf^  - ^(mf^  + 1 ) 

a(mf^  - a)(mf2  + 1)  ^ 


(5.37) 


and  is  plotted  in  the  a-3  plane  for  different  values  of 
boundary  ratios  e.|  and  £2*  we  notice  from  Figure 
5.8,  the  solid  boundaries  tend  to  destabilize  the  flow 
for  small  azimuthal  wave  numbers  as  they  are  brought  in. 
This  destabilization  is  at  first  sight  unexpected  since 
it  contradicts  what  would  be  expected  in  the  case  of 
continuous  profiles  as  the  walls  are  brought  sufficiently 
close  to  the  interface.  However,  the  mechanism  behind  it 
can  be  understood  by  considering  a two-fluid  system  between 

walls,  with  constant  density  and  velocity  in  each  region  in 
a two-dimensional  g ra v i ta i tona 1 flow  field  as  shown: 
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Figure  5.8  Effect  of  solid  walls  on  the  stability  boundary  for  ir 
subject  to  azimuthally  periodic  disturbances. 
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A well-known  result  for  interfacial  waves  in  the  absence 
of  surface  tension  for  given  values  of  wave  number  k,  ob- 
tained by  Greenhill  (1887),  is 

kp-|(U^  - c)^  cot  h(kh^)  + kp2(U2  - c)^  cot  h(kh2)  = g(p2  - P-;)  (5.38) 


Here  c is  complex  eigen  velocity  to  be  determined.  Con- 
sider long  wave  solutions  for  finite  wall  distances. 
Then  (5.38)  gives 


(U,  - c)^  (U2  - E)^ 

^1  hE|  ^2  h^  ^ 


or 


+ - 


a 


c = 


hih2 


(Ui  - 82)^  + (yp  + ^)9 


JL  + ^ 

hi  h2 


(5.39) 


is  the  density  ratio  at  the  interface.  The  destabiliza- 
tion of  moving  the  walls  inwards  can  be  observed  immedi- 
ately by  assuming  p.j  = p2  (i.e.,  a = 1)  and  is  similar 
to  what  we  observe  in  a rotating  case. 
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This  boundary  effect  for  the  present  rotating 
flow  will  diminish  as  m increases  and  eventually  it  will 
not  influence  the  stability  boundary  for  sufficiently 
short  azimuthal  waves. 

5.1.3  The  Non-rotating  Jets 

The  solution  of  two  concentric  non-rotating  jets 
with  different  densities  in  the  absence  of  solid  boundar- 
ies has  been  obtained  by  Alterman  (1961),  and  the  stability 
characteristics  for  some  cases  have  been  discussed.  We 
would  like  to  bring  up  here  some  of  the  interesting  re- 
sults concerning  the  bound  on  the  complex  eigen  velocity 
(c  = as  derived  in  the  discussion  of  the  general  cri- 
teria. 

For  = 9.^  = 0,  equation  (5.18)  can  be  arranged 
as 

E^(W^  - c)^  - aE2(W2  - c)^  = -F(l  - m^  - k^)  (5.40) 

where 
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Solving  for  the  complex  c,  one  finds 

- aE2W2  ±JaE^E^{Vi^  - - (E^  - oE2)r(i  - -x^) 


(5.41) 

The  condition  for  stability  is 

aE^E2(W2  - - (E^  - aE2)r(l  - - k^)  > 0 (5.42) 


Since 

m'  2'  — m'  ' - m'  r 

K'(S)  > K'(<)  > K'(k,) 
m2  — m — m l 


for  <2  — — '^1 


it  follows  that 

E^  > 0 E2  1 0 


We  notice  immediately  that  (5.42)  is  always  violated, 
i.e.,  instability  occurs  whenever 

r(l  - m^  - K^)  <0 


Surface  tension  has  stabilizing  effects  except  for  m = 0 
with  K<1,  and  the  flow  is  always  unstable  in  the  absence 
of  surface  tension. 


For  unbounded  flow. 


I (k) 
m ^ ' 

kI  ' (k) 
m 


K (k) 
m ~ ' 

K K ' ( K ) 
m 


equation  (5.40)  becomes 
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„ kK  ' ( k)  kI  ' ( k)  , , 

“<“2  - “l>  * - k“^)  < 0 (5.43) 

^m'  ^ m'  ' 

This  result  was  originally  due  to  Alterman  (1961). 

For  maximum  growth  rates,  we  differentiate  the 
imaginary  part  of  the  wave  speed  c for  unstable  solutions 
with  respect  to  the  density  ratio  a in  (5.41)  and  obtain 


1 'max 


2*^  1^1 


W, 


r(l 


li^l 


m^  - K^) 


(5.44) 


with  the  corresponding  real  wave  speed 


i[(W^  + W2)  + 


- 


E^(W^  - W^) 


and  density  ratio 

E,  [E,(W,  - - r(l  - - <^)] 

’ ^2  [E,(W,  - + r(l  - 


If  we  plot  the  unstable  solutions  in  the  complex  c plane, 
one  of  the  very  interesting  results  observed  is  that  the 
complex  eigen  values  lie  exactly  on  a certain  semi-circle 
governed  by  the  equation 


Where  the  positive  sign  corresponds  to  W.j  > and  the 
negative  sign  corresponds  to  W.|  < W2 . The  locus  of  maximum 


no 


growth  rate  is  found  to  be  a straight  line  in  the  complex  c 
plane  with  the  equation; 

Ci  - Cr  = - W^l  - (W^  + W^)]  (5.46) 

For  W-|  = W^,  the  only  instability  mechanism  is  due  to  the 
surface  tension  existing  at  the  interface,  the  flow  will 
only  be  unstable  when  m = 0 and  k < 1.  This  situation  was 
first  discovered  by  Rayleigh  (1892b)  for  unbounded  flows. 

Figure  5.9  shows  the  surface  tension  effect  on  the 
unstable  solutions  subject  to  axisymmetric  disturbances 
for  different  values  of  k in  the  complex  c plane.  A larger 
axial  wave  number  corresponds  to  a smaller  semi-circle. 

All  these  semi-circles  will  eventually  shrink  to  a point 
at  c^.  = 0 and  c^  = W2,  where  the  density  ratio  a is  equal 
to  infinity,  and  the  flow  is  stable.  If  we  look  back  to 
(5.42),  we  notice  that  the  damping  effects  caused  by  sur- 
face tenison  at  the  interface  will  be  significant  for  suf- 
ficiently large  axial  and  azimuthal  wave  numbers.  In  such 
cases,  the  curvature  of  the  disturbed  interface  will  be 
large  enough  to  produce  damping  forces  to  stabilize  the 
onset  of  instability.  Therefore,  the  unstable  solutions 
for  sufficiently  short  wave  lengths  will  not  exisit  due  to 
these  stabilizing  effects. 
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wave  speed  for  axisummetric  disturbances. 
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In  the  absence  of  the  surface  tension,  equation 
(5.45)  gives 

- ]-(W^  + W^)}^  + = {^(W^  - W^)}^  (5.47) 

This  is  also  very  interesting  as  the  above  equation  is 
exactly  the  semi-circle  bound  on  the  complex  c given  by 
(3.29).  All  the  solutions  (in  this  case,  they  are  all 
unstable)  lie  not  inside,  but  exactly  on,  a semi-circle 
with  the  range  of  the  axial  velocity  as  diameter.  A 
question  may  be  raised  at  this  point  how  the  exactness 
between  the  predicted  bound  and  the  solution  can  ever  be 
achieved  even  in  the  presence  of  the  density  and  boundary 
effects.  In  fact,  this  can  be  explained  by  looking  back 
to  the  derivation  of  the  semi-circle  bound  on  the  complex 
wave  velocity  for  non-rotating  jets. 

Combining  equations  (3.28),  one  obtains 

^2  ^2 

0 £ / (W  - a)(W  - b)Qdr  = {[c^  - ^(a  + b)]^  “ b)]^l  / Qdr 

R,  R, 

(5.48) 

For  the  flow  under  consideration,  the  axial  velocity  pro- 
file is 

W = + (W^  - W.J  )H(r  - R) 


(5.49) 
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with  the  lower  bound  of  the  velocity  W-j  = a and  the  upper 
bound  of  the  velocity  W2  = b,  assuming  W2>  . H(r  - R) 

stands  for  the  Heaviside  unit  step  function.  From  (5.48) 

R2 

! (W  - W, )(W  - W.jqdr 


P 

R-  c R+c  ^2 

= 1im(  / + / + / > (W  - W,)(W  - W2)Qdr 

e-»-o  R-£  R+e 


R+e 

= lim  / [{W2  - W^)H(r  - R)][(W^  - W2)  + {W2  - W^)H(r  - R)]  Qdr 

e->oR-e 


^R+e 

= lim(W2  - /[H(r  - R)  - 1]  H(r  - R)  Qdr  ^ 0 (5.50) 

£->■0  R-£ 

In  this  case,  the  inequality  (5.48)  becomes  an  equality, 
and, accordi ngl y , all  the  unstable  solutions  with  the  same 
axial  velocity  profile  must  lie  exactly  on  the  semi-circle 
described  by  (3.29)  or  (5.47)  irrespective  of  either  the 
density  distributions  or  the  boundary  effects.  The  solutions 
for  very  small  or  very  large  density  ratios  lie  on  the  two 
ends  of  the  semi-circle  representing  the  stabilizing  effect 
of  density  on  the  growth  rates. 
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Figure  5.10a  shows  the  maximum  growth  rates 
(c^.  = ^ IW2  - W.jl)  for  different  values  of  k and  m for 
unbounded  flows.  Solutions  for  axisymmetric  disturbances 
are  characteristically  different  from  those  for  non- 
axisymmetric  ones.  For  sufficiently  large  azimuthal  or 
axial  wave  numbers,  the  curves  for  maximum  growth  rate 
approach  a = 1,  i.e.,  the  flow  will  be  most  unstable  when 
the  densities  of  the  fluids  approach  each  other  at  the 
interface.  This  behavior  is  similar  to  that  of  the  Kelvin- 
Helmholtz  instability  encountered  in  a two-layer  hetero- 
geneous shear  flow  (Lamb  1932,  §232).  The  maximum  growth 
rate  for  sufficiently  large  wave  numbers,  or  in  the  absence 
of  the  gravitational  and  surface  tension  effects,  occurs 
when  the  densities  in  the  two  layers  are  equal. 

The  effects  of  the  solid  boundaries  for  axisym- 
metric disturbances  are  plotted  in  Figures  5.10b,c,d.  The 
inner  boundary  travelling  outwards  will  move  the  curve  for 
maximum  growth  rate  up  in  the  a - k plane  while  the  outer 
boundary  travelling  inwards  will  move  it  down.  If  we 
bring  in  both  boundaries  sufficiently  close,  the  curve  for 
maximum  growth  rate  will  approach  the  line  a = 1.  Indeed, 
the  flow  will  always  be  stable  for  all  azimuthal  wave  num- 


bers if  k = 0. 
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0 1 2 3 4 5 

K = kR 

Figure  5.10a  Curves  of  maximum  growth  rate  for  different 
values  of  wave  numbers. 


Figure  5.10b  Curves  of  maximum  growth  rate  for  different 
values  of  inner  boundary  ratio 


3 


Figure  5.10c  Curves  of  maximum  growth  rate  for  different 
values  of  outer  boundary  ratio  £2- 
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5 . 2 A Rotating  Jet  Core  with  Continuous  Density  Distri- 
bution Embedded  in  a Potential  Vortex  of  Different 
Density 

The  stability  of  a potential  vortex  with  axial 
motions  was  first  studied  analytically  by  Ponstein  (1959) 
including  the  surface  tension  and  viscosity  of  the  liquid, 
and  the  inertial  effects  of  the  ambient  air.  He  found 
that  in  some  cases  non-axi symmetr i c perturbations  grew 
more  rapidly  than  axisymmetric  ones  and  instabilities  took 
place  if  the  wave  numbers  in  both  directions  were  less 
than  some  critical  values.  Also,  the  influence  of  the 
inertial  effects  of  the  ambient  air  on  the  solid  jet  might 
be  neglected  if  the  density  of  the  air  was  small  compared 
to  that  of  the  liquid  and  provided  that  the  velocity  field 
was  continuous  across  the  interface.  Michalke  and  Timme 
(1967)  approximated  a single  vortex,  as  it  occurs  in  a 
free  boundary  layer,  by  an  axi symmetri  cal  vortex  model 
of  a homogeneous  fluid  subject  to  inviscid  stability  an- 
alysis. One  of  the  important  features  they  pointed  out 
is  that  the  two-dimensional  azimuthally  periodic  pertur- 
bations are  apparently  more  strongly  amplified  than  the 
corresponding  three-dimensional  ones.  This  has  been  ob- 
served experimentally  by  Weske  and  Rankin  (1963)  in  a 
special  arrangement.  However,  in  one  of  the  three  special 
types  of  vortex  they  examined,  they  failed  to  consider 
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the  centrifugal  effects  arising  from  the  discontinuity 
of  velocity  field  across  the  interface,  which  introduced 
error  into  their  analysis.  Later,  Uberoi  et  a 1-  ( 1 972  ) 
looked  at  a rotating  jet  embedded  in  a single  vortex  with 
different  densities  and  obtained  approximate  expressions 
for  the  complex  wave  numbers  valid  for  very  long  and  very 
short  axial  wave  lengths.  Recently,  Lessen  et  a 1.  ( 1 973  ) 
investigated  a similar  problem  in  the  limiting  case  of 
homogeneous  fluid  and  procured  numerical  solutions  for  a 
non-rotating  and  rotating  jet  core  with  no  discontinuity 
in  azimuthal  velocity  across  the  interface. 

In  the  following  example  we  investigate  the  sta- 
bility of  a cylindrical  potential  vortex  surrounding  a 
rotating  jet  core  with  continuous  density  variation. 

The  flow  under  consideration  consists  of  two  re- 
gions of  fluid  with  an  interface  at  r = R and  the  inner 
and  outer  solid  boundaries  at  R^  and  R2  respectively. 

The  velocity  and  density  distributions  in  regionsl  and  2 
are  given  as 

' 

< W = W-|  for  R-j  £ r < R 

^0  = PK^)"" 


- W = W2 
pQ  = P2 


1 1 9 


for  R < r £ R2 


a„(r)  P„(r) 


Figure  5.11  Velocity  and  density  distributions  for 
a rotating  jet  core  embedded  in  a 
potential  vortex. 

Here  ^2’  ^1’  ^2  ^ constants. 

Before  proceeding  to  the  analysis  of  the  flow,  we 
would  like  first  to  discuss  the  stability  characteristics 
concerning  the  velocity  and  density  profiles  in  region  1. 

Due  to  the  difficulty  in  solving  equation  (2.11) 
for  continuous  velocity  and  density  distributions,  analytical 
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solutions  have  up  to  now  been  confined  to  broken-line 
type  profiles  in  velocity  and  especially  step  function 
type  profiles  in  density.  The  smoothly  varying  density 
distribution  we  are  now  discussing  is  a break  from  these 
confines. 

Equations  (2.9)  and  (2.10)  for  the  profiles  in 
region  1 reduce  to 


n^Dp^  + ^p-j  = -if2^p^(^)  [(a  + 4)  - n^]u^ 
2 ^ 

(k^  + ^^)p^  = i^^,P-,(^)  (^  - n^D*)u^ 

r 


where 


n 


1 


+ m 


(jO 


(5.51 ) 


(5.52) 


Eliminating  u-j  from  the  above  two  equations,  we  obtain 


+ q.k^)Pi  = 0 


(5.53) 


where 


4 + 0 


V = 


2mo 
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a 

Upon  letting  = r^y,  equation  (5.53)  can  be  transformed 
into  a modified  Bessel  equation 

2 

D^Y  + -DY  - + q^k^)Y  = 0 (5.54) 

with  solutions  in  the  form  of  the  modified  Bessel  function 
of  the  first  and  second  kinds,  namely 

Y = A*I^(qkr)  + B*K^(qkr)  (5.55) 

where 


V 


+ oA  + (|)^ 


The  solution  for  the  perturbation  pressure  can  now  be 
written  as 

a 

p,  = -1n^p,(^)“(n^  - 4 - o){r  + Bk'"  K ( qkr- ) ] ) 

(5.56) 

where  k"^  and  k'^  are  normalizing  factors  in  order  to  allow 
k ^ 0 for  V / 0. 

Substituting  (5.56)  into  (5.51),  we  obtain  the  solution 
for  the  perturbation  velocity  in  the  following  form. 
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“l  = N,r 


- ' -v.2m  , o . 


„ qkrK' ( kqr) 

V r2m  , a . ^ V ' ^ ' 


(5.57) 


Here  we  consider  the  profile  in  region  1 to  be 
within  the  annulus  of  two  solid  walls  at  r = R-|  and  R2 
for  the  time  being.  Some  of  the  stability  conditions 
for  this  continuously  varying  density  profile  can  be  re- 
vealed by  making  use  of  the  boundary  conditions  that  the 
perturbation  velocity  must  vanish  at  r equal  to  R-j  and 
R2.  This  yields 


- qKiI'CqK,) 

[2S1  t o + 1,^,  \]1  (q<, ) 

•-n^  2 V^i^  ^ ^ 


[^  + I + 

"1  2 


V 2 


_ qK,K'(qK,) 

+ o ^ .1  V 1 1 ( ) 

•-n-,  2 K^Iqi^  1 


„ qK„K'(qiCo) 

r2m  , a ^ 2 2 -,y  , \ 

[h^^2  * R;;Tql7p'-'%'‘’^2' 


(5.58) 


Since  the  amplification  factor  w is  inside  not 
only  the  argument,  but  also  the  order  of  the  modified 
Bessel  functions,  analytical  solutions  are  not  possible 
except  for  special  conditions.  Two  particular  cases  will 
be  discussed  and  compared  with  the  stability  criteria: 

(a)  For  axisymmetric  disturbances,  equation  (5.58) 
can  further  be  reduced  to 
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Ij,  (q<2)  (qicj) 


(5.59) 


where 


u = I 1 - al 

An  asymptotic  expansion  of  the  Bessel  functions  in  (5.59) 
yields 


Recalling  q in  (5.53)  and  solving  for  w , one  finds. 


0)  will  be  complex  if  o < -4,  which  is  exactly  what  the 
Rayl ei gh-Synge  condition  predicts.  In  fact,  the  axi- 
symmetric  modes  will  be  stable  for  all  axial  wave  lengths 
whenever  a ^ -4,  according  to  that  condition. 

(b)  For  small  axial  wave  numbers,  the  modified  Bessel 
functions  take  their  algebraic  forms  and  (5.58)  is  re- 
ducible to  a very  simple  form 


g2q(<2  " ^1 ) 


e 


2aTT  i 


a = 0,  1 , 2 . . 


(5.60) 


a 


=0,1,2.  . 


R.  2v 


1 24 

2aiTi 


The  above  equation  is  equivalent  to  a quadratic  equation 
in  n.| 


+ 


2mani 


2 


tn  o 


0 


or 


0) 


, *^1  . ma  ± 

\ " — 


^/o  { 0 + a”} 


A 


(5.61) 


where 


A 


log(^) 


Since  the  term  inside  the  curly  bracket  is  positive 
definite,  the  stability  condition  depends  only  on  the 
sign  of  a,  the  density  characteristic  parameter.  Equa- 
tion (5.61)  is  exact  for  k - o, saying  that  negative  a 
(representing  negative  density  gradi  ent)  i mpl  i es  instability, 
while  positive  a (representi ng  positive  density  gradient) 
implies  stability,  which  is  exactly  as  predicted  in  Section 
(3.3.2),  Chapter  3.  The  corresponding  A in  thi-s  case  is 
exactly  equal  to  1/4.  For  this  reason,  it  is  concluded 
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that  the  upper  bound  on  A for  which  instabilities  are 
possible  to  exist  can  not  be  improved  anymore.  Besides, 
it  is  well  known  that  the  solid  body  rotation  for  homo- 
geneous fluids  (corresponding  to  a = 0)  is  sta.ble  against 
arbitrary  disturbances.  The  above  result  clearly  reveals 
that  the  two-dimensional  azimuthal  modes  are  the  most  un- 
stables ones,  and  the  difference  between  the  onset  of  in- 
stability for  axisymmetric  modes  and  non-axi symmetri c ones 
is  quite  significant. 

With  the  understanding  of  the  stability  behavior 
of  the  flow  in  region  1,  we  now  proceed  to  the  analysis  of 
the  two-region  flow  under  consideration. 

Consider  only  positive  o for  centrifugally  stable 
profile  in  region  1.  The  corresponding  solutions  for  per- 
turbation pressure  and  velocity  are  previously  given  in 
(5.56)  and  (5.57). 

For  the  profiles  in  region  2,  equation  (2.11)  can 
be  reduced  to 


D[-o — — y-^(Du„  + 


0 


(5.62) 


Let 


U2  = D(j) 


Integrating  (5.62)  and  adjusting  the  integration  constant. 
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one  finds 


D^(|) 


2 9 

(\  + k^)(j)  = 0 


<)>  = A2ln^(kr)  + B*K^(kr) 
or 

ijkr)  t B2k'”*'>K^(kr)] 


(5.63) 


(5.64) 


The  perturbation 
stituting  (5.64) 

2 ^2 
P2 


pressure  in 
into  (2.9) 


+ m(^)  - 


region  2 is  obtained  by  sub- 
and  has  the  form 


][A2k-"'l^(kr)  H-  B2k™K^(kr)] 

(5.65) 


Since  the  boundary  conditions  require  the  perturbation  ve- 
locity to  van  i sh  at  and  , the  solutions  in  region  1 
and  region  2 become 


^ n qkr I ' ( kqr ) 

= A.k-'^r  2 — I3i^]l„(l<qr) 


N,  "1 


"I  2 


V 


- 9m  ^ qkrK'(kqr) 

H [Ijn  t I + , V 

' "l  2 


N V » y 
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^ (^)  (n^  - 4 - a){r  ^[I^(kqr)  - H^K^(kqr)]} 


B2k'"[kK^(kr)  - H2kr(kr)] 


(5.67) 


Pz  - 


Hglmlk'-)] 


where 


„ qK, I ' ( qK, ) 

[2®  + o + 1..]I  (qK.) 

'-n^  2 I^(qK-|  ) 1 

r2m  , a , ) 


^;<*=2> 


The  secular  relation  for  stability  characteristics  nov; 
obtained,  by  matching  both  kinematic  and  dynamic  conditions 
at  the  interface,  is 


(n^  - 4 - o)[I',(qK)  - H,Kyq,c)] 


q kI ' ( q k)  ^ n 

f * - "itff  ^ 


1 


qKk’  (qK) 


2 2 


H2lm(K) 


"2  CTk7-  Ho  KlirKTy 


a 6 


1 - r,  (1  - m^  - K^) 


m 


m 


(5.68) 
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a,  6,  K and  have  the  same  definitions  as  those  in 
(5.18) . 

/s  /N 

In  the  absence  of  boundaries,  both  H-|  and 
will  go  to  zero  and  (5.68)  reduces  to 

7 7 ^ ^ 2 

"2  ■ 

(1  - m^  - K^)  (5 

Uberoi's  result  can  be  obtained  by  letting  a = 0 and 
in  the  above  equation. 

Several  particular  cases  of  secular  relation 
will  now  be  investigated; 


5.2.1  Axi symmetri c Perturbations  (m  = 0) 


Since 

i;(z)  = 

‘v  - - 

k;(z)  = 

-Kv  - 

- 

equation  (5.68) 

f or  m = 

0 can  be  written  as 

''7  ''2 

qF,N^  = 

7 2 

- 1 - l'i(l 

- r^)] 

'1 


2 I^{qKT-‘ 


r2m 

"1 


.69) 

= 0 
(5.68) 


.70) 


where 


1 29 


i^(qK) 

■v  - 


1 (q^i ) 


VTTT^ 

(qK-, ) 


K^(qK) 


kTTTT^  T^q^T 


'V  - 1 


Iq^) 


(qK) 


K^(k) 

K^{kJ 


1 - 


K-|(k2) 

171^ 

K 1 { K 2 ) 17*^^ 

I7K2)  171^ 


If  there  is  no  axial  velocity  difference,  i.e.,  W-j  - W2> 
equation  (5,70)  has  a simple  form 


(c 


[a3^ 


1 - r^i 


K^)] 


k ( q F ^ + a F 2 ) 


(5.71  ) 


Since  F2  and  the  real  part  of  7 are  positive  definite  for 
Ki  < K £ the  stability  condition  is  now  determined  by 

the  terms  inside  the  square  bracket,  and  the  stability 


boundary  for  axisymmetric  disturbances  and  zero  axial 
flow  gradient  is  again 


- 1 = r^l  - K^) 


(5.72) 


This  relation  for  zero  surface  tension  is  also  obtainable 
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by  integrating  the  Rayl ei gh-Synge  condition  across  the 
interface. 

5,2.2  Azimuthally  Periodic  Perturbations  (k  = 0) 

For  0=0  and  small  axial  wave  numbers,  equation 
(5.68)  can  be  written  in  an  algebraic  form  by  the  small 
argument  expansion  of  the  modified  Bessel  functions.  The 
secular  relation  is  reducible  to  the  quadratic  equation 


(5.73) 


with  f-j  and  f^  defined  in  (5.29). 


Equation  (5.73)  allows  one  to  solve  for  w and 


yields 


[f^(kll^  + m^l^)  - + af2(kW2  + -JhjH 


f-|  + af2 


with 


A = -afif2[Y-|  + - 1)  + ctf2)^f 


+ 


r,(i  - m^  - K^)]} 


The  stability  criterion  corresponds  to 
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A ==  0 (5.75) 

Several  phenomena  obtained  for  special  cases  follow: 

( a ) Effects  of  axial  wave  numbers 
For  unbounded  flow,  the  stability  boundary  with- 
out surface  tension  for  k = 0 is 

a[m(6  - 1)  + 1]  - (1  + ct)[1  + m(a3^  - 1)]  = 0 (5.76) 

Figure  5.12  shows  the  stability  boundaries  for  different 
values  of  m for  k = o.  These  curves  do  not  intersect 
each  other  and  those  with  larger  wave  numbers  correspond 
to  the  less  stable  modes.  Such  azimuthal  modes  will  am- 
plify in  time  and  the  stable  region  will  eventally  shrink 
into  the  line  B = 1>  ot  ^ 1.  The  lines  aB  = 1 correspond- 
ing to  stability  boundaries  for  axisymmetric  modes  are 
plotted  and  compared  to  those  for  azimuthally  periodic 
ones.  Once  again,  it  is  observed  that  those  stability 
boundaries  lie  partially  in  the  region  stable  to  axisym- 
metric modes  in  the  a - B plane.  This  means  the  present 
flow  can  be  quite  stable  for  axisymmetric  perturbations 
yet  unstable  to  non-a xi symmetr i c ones.  These  results  also 
show  that  the  Rayl ei gh-Synge  criterion,  correspondi ng  to 
aB^  - 1 in  the  present  example,  is  not  valid  when  non- 
axisymmetric  perturbations  are  taken  into  account,  and 
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Figure  5.12  Stability  boundaries  for  axis  mmetric  disturbances  and 

azimuthally  periodic  disturbances  of  different  values  of 
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that  two  dimensional  m modes  are  most  unstable. 

( b ) Effects  of  axial  flow  difference  for  small  k 

The  influence  of  large  axial  flow  difference  (i.e., 

IW2  " W-|l  large)  is  still  of  importance  even  for  small  axial 
wave  numbers.  Equation  (5.75)  without  surface  tension  and 
solid  boundaries  takes  the  form 

2 ? 

a[Yi  + m(B  - 1)  +1]  - (1  + ot ) [ 1 + m(a6  - 1)]  = 0 (5.77) 

The  stability  boundaries  are  plotted  in  Figure 
5.13.  The  behavior  of  those  curves  under  the  influence 
of  axial  flow  difference  resembles  that  of  the  two- region 
flow  with  constant  angular  velocities.  Positive  y-j  has 
destabilizing  effects  when  B is  positive  and  stabilizing 
effects  when  B is  negative;  and  vice  versa  for  negative 
■Y.|  . Figure  5.14  shows  the  stability  boundaries  for  dif- 
ferent values  of  m in  the  a - y-|  plane,  the  angular  veloc- 
ity being  continuous  across  the  interface.  All  the  stabil- 
ity boundaries  pass  through  the  point  a = 1 , y-|  = - 1 . 

The  decrease  of  density  and  the  increase  of  axial  flow 
difference  both  destabilize  the  flow. 

( c ) Effects  of  boundaries 

The  stability  boundary  described  by  equation  (5.75) 
in  the  absence  of  surface  tension  is 
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for  small  axial  wave  numbers. 
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Figure  5.14  Stability  boundaries  of  continuous  angular 
velocity  for  small  axial  wave  numbers. 
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+ m{6  - 1 )]  - (f^  + a + m(a3^  - 1 )]  = 0 

(5.78) 


Figure  5.15  reveals  that  the  flow  is  destabilized  by  the 
inward  moving  walls  for  m = 2.  This  destabilization, 
contradicting  physical  expectations  especially  when  the 
solid  boundaries  are  sufficiently  close,  is  a result  of 
the  assumed  discontinuity  arising  in  the  velocity  field 
as  explained  previously  in  the  two -region  flow  for  con- 
stant velocities.  This  phenomenon  diminishes  as  m in- 
creases and  the  stability  boundary  will  be  unaffected  by 
the  walls  for  sufficiently  short  azimuthal  waves. 


5.2.3  A Non-rotating  Jet  Core  with  Density  Variations 
Surrounded  by  a Homogeneous  Potential  Vortex 

For  = 0,  equation  (5.68)  yields 

2 ^ 2 p 2 2 

E^(kW^  - to)  - aE2(kW2  + (o)  = a^pCl  - Tp  ( 1 - m -k)] 

(5.79) 


where 


K I ' ( kT) 

[f  + 

I^(k) 


_ kK ' ( K ) 

]I^(k)  - H^[|  + 


^(r) 


137 


b 


OsJ 

II 


Figure  5.15  Effect  of  solid  walls  on  the  stability  boundary  for 
subject  to  azimuthally  periodic  disturbances. 
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Equation  (5.79)  is  a quadratic  equation  for  w and  can  be 
solved  as 


U) 


E.|(kWi)  - aE2(kW2  + 


E 1 ~ ot  E 2 


^2  ^-aE^  E2  (y2  " “(^1  ■ aE2)[1  " r2(l  - - k^)] 

^ ' - /V  /N 

E-j  - aE2 


with 


^2 


k(W2  - W^) 
^2 


(5.80) 
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The  flow  will  be  stable  whenever 

E^E2(y2  + + (E]  - o‘E2)[1  - r2(l  - - k^)]  > 0 (5.81) 

Since  > 0,  E2  1 0 (see  Appendix),  both  the  axial  veloc- 
ity difference  and  the  azimuthal  wave  numbers  have  de- 
stabilizing effects  provided  they  are  of  the  same  sign. 
Surface  tension  always  stabilizes  the  flow, except  when 
m = 0 and  k < 1 . 

For  maximum  growth  rates  of  the  unstable  solutions, 
one  differentiates  the  imaginary  part  of  (5.80)  with  re- 
spect to  a and  obtains 


(jO. 


ma  X 


[1  + 


1 - r^d  - m^  - K^) 

E2(Y2 


] (5.82a) 


with  the  corresponding  real  part  of  w 


(jj 


^[kW.|  + kW2  + mfl2  + ^2 


1 - T2(1  - m^  - K^) 

E2(Y2  + 


] (5.82b) 


and  density  ratio 


E,  E (yo  + m)  + 1 
. -(^)  2 2 


T2(1  - m^  - k2) 


E2  E2(y2  + 


2 2^ 
m - K ) 


a 


1 + ^2^^ 


(5.82c) 
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It  is  interesting  to  note  that  all  the  unstable  solutions 
for  the  flow  under  consideration,  irrespective  of  the 
types  of  perturbations,  the  boundary  effects  or  the  sur- 
face tension  existing  at  the  interface,  lie  exactly  on 
a certain  semi-circle  governed  by  the  equation 

2 2 

{(jj  )}  +u?  = la).|  (5.83) 

' max  max 

The  positive  or  negative  sign  corresponds  to  the  positive 
or  negative  value  of  ^ respectively.  The  locus 

of  maximum  growths  lies  exactly  on  a straight  line  in  the 
compl  ex  0)  plane. 

We  will  now  discuss  some  of  the  stability  charac- 
teristics concerning  unbounded  and  bounded  flows. 


( a ) Unbounded  flows 


For  unbounded  flows,  both  H.|  and  in  (5.79)  go 


to  zero,  so  that  the  equation  yields 


(kW^  - oj)^ 


2 2v  n 

m - K ) J 


(5.84) 


This  flow,  in  the  absence  of  surface  tension,  density  vari- 
ations and  axial  velocities,  was  originally  investigated 
by  Michalke  and  Timme  (1967)  for  arbitrary  disturbances. 
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They,  however,  failed  to  consider  the  centrifugal  effect 
due  to  the  angular  velocity  jump  at  the  interface. 

In  this  particular  case,  the  neutral  stability 
condition  as  deduced  from  (5.81)  is 


2 K I ' ( K ) 

(m+  Y2)  - 


kK' (k) 
m ' ' 

K (k) 
m 


} = 


0 


and  stability  can  be  guaranteed  whenever 


a > 


2 

(m  + Y2)  + 


m 


i'i<) 

17^ 


(5.85) 


Figures  5.16  a,  b,  c show  the  stability  boundaries  for 
different  values  of  a and  m.  Both  the  increase  of  the 
relative  axial  velocity  and  the  decrease  of  the  density 
ratio  at  the  interface  have  destabilizing  effects.  The 
unstable  solutions  for  a = 0 and  m = 0 are  plotted  in 
Figure  5.16d,  which  are  all  exactly  on  a semi-circle 
governed  by  the  equation 

{c  - (W  + [c  . I )}  + c - Ic ■ I (5.86) 

' max  max 


where 
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Y2  = k{W2  - W^)/fl2 


Figure  5.16a  Stability  boundaries  for  m = 2 and  a - 0 
(symmetric  w.r.t.  Y2  " "2)- 


0 1 2 3 4 5 


Y2  = k(W2  - W^)/J^2 

Figure  5.16b  Stability  boundaries  for  m = 0 and  a = 2 
(symmetric  w.r.t.  Y2  ~ 0)- 
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Figure  5.16d  Effect  of  axial  wave  numbers  on  the  complex 
wave  speed  for  m = 0 and  a = 0. 
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with  the  locus  of  maximum  growth  rate  on  astraight  line 

c-c.=W,  (5.87) 

r 1 1 

Another  result  which  should  be  noted  is  that  all 
the  unstable  solutions  for  the  axisymmetric  case  lie  in- 
side a semi-circle  with  the  range  of  the  axial  velocity 
as  diameter.  It  gives  a solid  support  to  the  semi-circle 
theorem  derived  in  equation  (3.31). 

Perhaps  it  should  be  emphasized  once  more  that 
the  Rayl ei gh-Synge  criterion  is  no  longer  valid,  even 
only  for  axisymmetric  disturbances,  once  the  axial  veloc- 
ity gradients  are  taken  into  account.  Instabilities  still 
take  place  through  the  shear  layer  arising  from  the  axial 
flow  difference  at  the  interface. 

( b ) Bounded  flows 

For  simplicity  the  stability  boundary  described 
by  equation  (5.81)  for  bounded  flows  can  be  approximated 
by  the  lower  and  upper  limits  of  axial  wave  numbers 
(though  there  is  no  difficulty  in  evaluating  the  equation 
for  arbi trary  k) . 

For  k = 0,  the  stability  boundary  is  governed  by 


the  equation 


a = ( 1 - -~)  (f , V - j) 

m?2 


where 


1 + 


1 - 


e 


e 


2v 

1 

2v 

1 


(5.88) 


/s  I • c-  2 

^2  ” 1 - e2™ 

Stabilities  can  be  quaranteed  when  the  density  ratio  at 
the  interface  exceeds  that  of  (5.88)  for  a given  value  of 
parameter  a. 

Figure  5.17  shows  the  behavior  of  the  stability 
boundaries  subject  to  the  change  of  boundary  distances, 
both  the  density  ratio  a and  the  density  parameter  a serve 
as  a stabilizing  factor  to  the  flow.  The  destabilizing 
effect  of  moving  the  boundaries  inwards  increases  with 
descreasing  azimuthal  wave  numbers, and  occurs  when  the 
distance  between  two  walls  is  sufficiently  close. 

For  large  axial  wave  numbers,  the  stability  bound- 
ary of  (5.81)  can  be  written,  after  some  algebraic  manip- 
ulation , as 
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Figure  5.17  Effects  of  boundary  ratios  on  stability  boundaries 
for  azimuthal 1y  periodic  modes. 
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a = 


(5.89) 


where 


f 


* 1 t 

1 ■ 1 _ e-2(K  - •'!> 


f 


* 


2 


The  stability  boundaries  are  plotted  in  Figure  5.18.  The 
increase  of  density  gradient  stabilizes  the  flow.  The 
destabilization  effect  of  the  solid  walls  in  the  present 
case  is  not  that  significant  when  compared  with  that  for 
k = 0,  and  will  not  be  noticed  until  the  distance  between 
walls  is  quite  close. 

5.2.4.  Non-rotating  Jets  with  Radial  Density  Variations 


Equation  (5.68)  for  - 0 yields 


E^(W^  - c)^  - aE2(W2  - c)^  = -r(l  - m^  - k^) 


(5.90) 


R eca 1 1 that 
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Figure  5.18  Effects  of  boundary  ratios  on  stability 
boundaries  for  k = 5. 
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r = 


and  and  E2  are  previously  defined  in  (5.79) 

For  a = 0,  v=  m , equation  (5.90)  reduces  exactly 
to  (5.40),  i.e.,  two  homogeneous  concentric  jets  with  dif- 
ferent densities.  Since  both  (5.40)  and  (5.90)  have  ex- 
actly the  same  form  except  with  a slightly  different  def- 

/\ 

inition  between  E-j  and  E-j  which  are  both  constants,  all 
the  analysis  and  results  we  have  carried  earlier  can  of 
course  be  applied  here.  Therefore,  it  is  concluded  that 
the  stability  condition  for  the  flow  is 


aE^E^O^^  - W^)^  - (E^  - aE2)r(l  - m^  - kc^)  > 0 (5.91) 

which  is  always  violated  in  the  absence  of  surface  tension, 
despite  the  solid  boundary  effects.  In  addition,  all  the 
unstable  solutions  in  the  complex  c plane  lie  exactly  on 
a certain  semi-circle  satisfying  the  equation 

2 2 2 

{c  - (W„  + lc.|  )}  + {c.}  = lc-1  (5.92) 

max  max 

with  the  maximum  growth  rate  lying  on  the  straight  line 


"r  ' ■ “2I  ■ 
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where 


o[  I Wi 


ma  X 


r, ( 1 - - K^) 

^2  I — : ■] 


El  |Wi 


W, 


with  the  positive  or  minus  sign  depending  on  the  sign  of 
{W^  - W2}  respectively. 

In  the  absence  of  surface  tension,  all  the  solu- 
tions, which  are  always  unstable,  lie  exactly  on  a semi- 
circle in  the  complex  c plane,  with  the  range  of  axial 
velocity  as  diameter. 

The  density  for  maximum  growth  rates  (C:j  = ^ 

1 i es  on  the  line 


H,  - W2 


a = — 


El  [Ei(W^  - W2)^  - r^(l  - m^  - K^)] 
1 - (1  - m^  - K^)] 


[Ei(W 


(5.93) 


In  the  case  of  zero  surface  tension,  the  above  expression 
can  be  approximated  by  using  expressions  valid  for  small 
and  large  axial  wave  lengths,  i.e.. 


a = 


f,  M * (|)^ 


a 

2 


for  K <<  1 


(5.94a) 


mf , 


a = 


(^) 


^ - K f 

7 '^r  1 


(f)2 


+ ^(k 


★ 

K,)f, 


KK- 


for  tc  >>  1 


(5 . 94b) 
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where  f-j  , f 2 , f-|  and  f2  are  as  defined  in  (5.88)  and  (5.89) 
Figures  5.19a,b  show  the  curves  of  maximum  growth  rate  (c^. 

- W2 I ) subject  to  the  change  of  solid  boundaries  and 
the  density  parameter  a.  For  sufficiently  long  or  short 
axial  wavelengths  with  a = 0,  the  curves  for  maximum  growth 
rate  will  approach  the  line  a = 1. 


Cl  £2 

Figure  5.19a  Effects  of  boundary  ratios  on  the  curves  of 
PBximum  growth  rate  for  m = 0 and  small  k. 
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Figure  5.19b  Effects  of  boundary  ratios  on  the  curves  of 
maximuin  growth  rate  for  k = 5. 
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5 . 3 Rigid  Body  Rotation  --  Th ree- req i on  Solution 

Rigid  body  rotation  is  of  practical  interest  in 
the  newly  developed  counterflow  centrifuge  method  for 
separation  of  uranium  isotopes,  and  also  is  one  of  two 
velocity  profiles  for  which  analytical  solutions  in  terms 
of  well-known  functions  are  possible. 

Consider  three  regions  of  fluids  with  the  inter- 
faces at  R and  R + 6 in  a constant  angular  and  axial  veloc- 
ity field  [0,  W^].  The  density  distribution  in  each 

region  is  shown  in  Figure  5.20. 
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Figure  5.20  Density  distribution  for  rigid  body  rotation. 


Due  to  the  unboundedness  of  the  modified  Bessel  function 
of  the  first  kind  at  infinity  and  the  second  kind  at  the 
origin,  the  solutions  for  the  perturbation  velocity  and 
pressure  in  each  individual  region  found  in  Sectionsl  and 
2 in  this  chapter  are  simply 

\ 1 % ■■ 

,p^  = - 4)[A^k‘"’l|^(qkr)] 


qkri ' (qkr) 


0 0 V 


qkrK' (qkr) 
2 ■ — 


P2  =-lP2(^)”«o<"o  - '*  - + BjkPK^fqkr)] 


83^  [ K„(qk>T  ^ 


•'P3“^"o  - '>)[B3k"'K^(qkr)] 


(5.95) 


where 
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=7’ -4 


1 - 


4 + a 


W 


- = + m 


0) 

Q. 


After  matching  the  kinematic  and  dynamic  boundary  conditions 
at  both  the  interfaces,  we  obtain  the  secular  relation  for 
stability  in  the  absence  of  surface  tension. 


.2m  0 ^ 


I 
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D^]I^(qK) 


„ qKK'C^ 

r2m  o ^ V ^ 

L„  ^9 


D^]K^(qK) 


„ qhl'(qh) 

r2m  , a . ^ ^ 

*-  n 2 r(qh) 


U ^ 2 7(  q 'h  T 


D2]I^(qh) 
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(5.96) 


where 


D,  = 


a. 


Hq  - (a  + 4) 
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R 


1 56 


The  density  ratios  at  the  interfaces  are  defined  as 


“1 


P2 

Pi 


a2 


P20  + 

P3 


(5.97) 


From  earlier  discussions  we  know  that  the  azi- 
muthally  periodic  modes  play  an  important  role  in  flow 
stability.  This  suggests  that  we  first  look  at  the  case 
of  finite  m and  small  axial  wave  numbers. 


5.3.1  Azimuthally  Periodic  Disturbances  (k  = 0) 
Equation  (5.96)  for  k = 0 yields 


(a,  - 1)A+  m + aAj  - v)  (a,  - 1)A  + m + a,  (|  + v)  Ov 

L_ ! = = — ( 1 + c ) 

(a2  - 1 )A  - m + - v)  (a2  - 1)A  - m + + v) 

(5.98) 

where  A has  been  defined  in  (3.37a). 

Two  special  cases  of  density  profiles  are  to  be 
discussed  as  follows: 


( a ) Step  function  type  profiles 

For  o = 0,  V = m,  equation  (5.98)  can  be  rearranged 


f[(p3  - P2)(p2 


)]A^  + 


2mf)2(p3 


- P-|)A 


m^[f(p3 


f*2 ) ( P2 


- Pi 


or 


2p2(pi  + P3)]  - 0 
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A = m 


p2(p-|  "P3)iA/^(p3”P2^  '^3^^  "4fPi  P2(P3"P2)  (P2“P]  ^ 

f(p3  - P2)(P2  ■ P]) 

(5.99) 

where  oK 


f = 1 - (1  + e) 


-2m 


1 > f > 0 


If  (P3  - P2)(p2  - P])  < 0’  density  distribu- 
tion either  behaves  like  Figure  5.21a  or  Figure  5.21b. 
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Figure  5.21  Step-type  density  distribution. 


Recall  equation  (3.39) 


0) 

mJ2, 


A - 1 ± a/I  - a 


(5.100) 
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Since  we  have  proved  in  Section  (3.3.2)  that  A has  to  be 
real  in  the  Sturm-Liouvil  1 e system,  instabilities  will  be 
expected  whenever 

A > 1 (5.101) 


Combining  (5.99)  and  (5.101),  and  considering  only  the 
minus  sign  in  front  of  the  radical  in  (5.99)  for  physical 
reasons,  one  finds  that  the  condition  for  instability  is 

(m^  - l)[f(p3  - P2)(P2  - - 2mfp2(P3  - P2)(P2  ' " ''^Pl 

(m  + DP3]  >0  (5.102) 


which  is  always  true  for  the  density  distributions  presently 
considered',  therefore  they  are  always  unstable.  This  sit- 
uation is  expected  if  we  look  at  the  general  criterion  for 
rigid  body  rotation  subject  to  azimuthally  periodic  per- 
turbations. Both  density  profiles  in  Figures  5.21a,b  have 
a negative  density  gradient  inside  the  flow  region  and  in- 
stabilities would  be  expected  in  both  cases.  If  p-|  = 03* 

equations  (5.99)  and  (5.100)  can  be  solved  for  w with 
P2 


U) 


1 ±J]  - mf 

' '■  "/V — ' ■ 
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m 


(5.103) 
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where 


a 


(1  + a)  + 


4 a" 


2m 


[(1  + s)""'  - 1] 


Since  f is  always  greater  than  one  for  positive  values  of 
m,  the  flow  is  always  unstable  except  for  a = 1,  i.e., 
rigid  body  rotation  for  homogeneous  fluid.  The  growth 
rate  versus  the  density  ratio  between  the  middle  ring  of 
fluid  and  the  surrounding  one  is  plotted  in  Figure  5.22. 

The  growth  rates  increase  in  both  directions  when  the  den- 
sity ratio  departs  from  the  line  a = 1 which  is  the  only 
stable  region  for  all  modes  in  the  - aplane.  The  higher 
azimuthal  wave  numbers  correspond  to  the  less  stable  modes; 
the  narrower  the  middle  ring  of  fluid  becomes  the  lower  the 
growth  rates  are.  When  the  density  of  the  fluid  in  the 
middle  region  is  sufficiently  large  or  small  in  relation 
to  that  of  the  surrounding  one,  the  flow  will  be  most  un- 
stable and  possess  very  large  growth  rates. 

If  (P3  - P2)(P2  - P])  = 0 (either  p-,  = P2  or 
P2  ""  p3),the  density  distribution  reduces  to  a two-region 
type  of  profile,  and  equations  (5.99)  and  (5.100)  can  be 
sol ved  for  w , i.e.. 


0) 


= m - 


1 - a + ^(a  - 1 ) [ ( m + 1 )a  + (m  - 1 ) J 


(5.104) 
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Unstable  modes  will  take  place  whenever  a < 1.  This  re- 
sult is  also  deducible  from  (5.33)  for  6=1.  Note  that 
the  unstable  solutions  for  the  complex  w in  (5.103)  and 

(5.104)  lie  exactly  on  a semi-circle  governed  by  the  equa- 
tion 


This  instability  characteristic  has  been  previously  dis- 
cussed in  the  general  criteria  in  Section  (3.3.2)  for 
rigid  body  rotation  subject  to  angul  arl  y per i odi c pertur- 
bations. 


If  a.|  = = 1>  the  density  will  be  continuous 

across  the  interface.  Equation  (5.98)  reduces  to 


If  the  middle  region  is  narrow  enough,  such  that  e <<  1 
and  (1  + e)^^=1  + 2ve,  equation  (5.105)  yields 


( b)  Continuous  density  distributions 


(5.105) 


0 


(5.106) 


For  positive  values  of  o,  the  density  profile  is  an 
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increasing  function  of  radius  and  the  flow  is  stable. 

For  negative  values  of  o,  eguation  (5.106)  reveals  that 
the  wave  velocity  will  be  complex,  and  hence  that  the 
flow  is  unstable.  Figure  5.23  shows  the  complex  wave 
speed  versus  the  density  factor  o for  e - 0.1.  The 
smaller  a is,  the  higher  the  growth  rate,  and  the  slower 
the  unstable  waves  will  propagate.  The  complex  wave 
velocity  lies  exactly  on  a semi-circle  as  forecast  in  the 
general  criteria  for  rigid  body  rotation  subject  to  azi- 
muthally  periodic  perturbations.  For  finite  values  of 
a,  the  solutions  for  the  complex  w lie  exactly  on  the 
right-=harld  side  of  the  semi-circle  as  shown  in  Figure  5.23. 


5.3,2  Arbitrary  Disturbances  for  Constant  Density  in  the 
Middle  Region 

Equation  (5.96)  for  a = 0 can  be  written  as 
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(5.107) 
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Figure  5.23  Growth  rates  for  azimuthally  periodic  pertur- 
bations (c  = 0.1). 
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wi  th 


The  stability  of  this  three-region  flow  for  rigid  body 
rotation  in  the  lower  limit  of  the  axial  wave  numbers 
has  just  been  discussed.  In  the  following  we  will  in- 
vestigate the  stability  behavior  of  the  present  flow  for 


large  axial  wave  numbers. 

Equation  (5.107)  for  short  axial  waves  can  be  ex- 
panded according  to  the  asymptotic  forms  of  the  modified 
Bessel  functions  to  yield 


{nf  - [2  + 


„ 1 - a,  2 2 


1 


[2  ij''*"  (t^>  “ 

(5.108) 


The  stability  condition  depends  on  the  density  ratio  at 
the  corresponding  i nter face, and  the  stability  boundaries 
are  independent  of  the  azimuthal  wave  numbers,  since  only 
the  unstable  modes  for  short  axial  wave  lengths  are  under 
consideration.  In  fact,  if  the  wave  lengths  are  short 
compared  to  the  thickness  of  the  middle  ring  of  the  fluid, 
instability  phenomena  will  be  confined  to  the  neighborhood 
of  the  interfaces  and  the  flow  will  not  be  affected  when 
it  is  sufficiently  far  away  from  the  interfaces.  Figure 
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5.24  shows  the  growth  rate  as  a function  of  density  ratio 
for  a-^  = = a,  i.e.,  p-j  = P2-  Since  instabilities  take 

place  at  the  interface  r = R only  when  a < 1 and  at  the 
interface  r = R + 6 only  when  a > 1,  the  conditions  govern- 
ing instability  following  from  (5.108)  are 


V4  + - 2 

V 4 + [k(1  + - 2 


for  a < 1 


(5.109) 


for  a > 1 


The  behavior  of  the  growth  rates  in  this  case  resembles 
that  of  azimuthally  periodic  ones.  The  growth  rates  for 
a > 1 increase  with  increasing  axial  wave  numbers  as  well 
as  the  increasing  width  of  the  middle  ring  of  the  fluid. 
The  unbalanced  centrifugal  pressure  produced  by  the  middle 
region  directly  affects  the  onset  of  instabilities.  The 
growth  rates  for  a<  1 depend  only  on  the  centrifugal 
pressure  balance  at  the  interface  r = R,  since  only  short 
wave  lengths  are  considered.  These  unstable  modes  will 
exist  for  all  values  of  a except  along  the  line  a = 1, 
i.e.,  solid  body  rotation  for  homogeneous  fluid,  which  is 
stable  against  any  disturbances. 
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ct  ^2!  P] 

Figure  5.24  Growth  rates  for  large  axial  wave  numbers 
(independent  of  m) 
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5 . 4 Non-rotating  Jets  --  Three-region  Solution 

We  know  that  the  instability  mechanism  for  non- 
rotating jets  is  a purely  shear  effect  between  two  ad- 
jacent fluid  layers,  and  that  the  particle  exchange  in- 
volves no  potential  energy  change  since  the  gravitational 
effects  of  the  fluid  are  neglected.  The  sufficient  con- 
dition for  stability  in  (3.19b)  fails  to  yield  any  infor- 
mation because  of  the  absence  of  stabilizing  effects 
existing  in  the  flow  field.  On  the  other  hand,  the  govern 
ing  stability  equation  for  non-rotating  jets  resembles 
Rayleigh's  equation  and  yields  only  a weak  necessary  con- 
dition for  instability.  As  a result  of  these  difficulties 
neither  stabilities  nor  instabilities  can  be  predicted  in 
advance.  However,  unstable  solutions  would  be  expected 
if  discontinuities  exist  in  the  axial  velocity  field. 

The  kinetic  energy  carried  by  the  fluid  particles  in  the 
neighborhood  of  the  discontinuous  sheet  will  suddenly  have 
a sharp  change  and  easily  achieve  the  particle  exchange 
due  to  the  absence  of  stabilizing  factors.  We  now  con- 
struct two  examples  to  demonstrate  the  instability  char- 
acteristics of  the  non-rotating  jet  flows. 

5.4.1  Three  Concentric  Constant  Axial  Jets 


The  first  example  for  non-rotating  jets  we  would 
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like  to  investigate  is  a three-region  step-type  velocity 
and  density  profile  with  the  interfaces  at  r = R and  r = 

R + 6 , i . e . , 


rW 


1 


r^l  for  0 < r < R 


W(r) 


P„('")  = P 


for  R < r < R + 6 
for  R + 6 < r < °° 


W(r) 


Figure  5.25  Step-type  distributions  for  three  concentric 
jets . 


Equation  (2.11)  for  zero  azimuthal  velocities  yields 
^ „ D*(p„EDW) 

D(p^ED  u)  - [p^k^  + — ^ = 0 (5.110) 
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For  constant  axial  velocities  and  densities,  equation 
(5.110)  becomes 


2 * 

D[(  r..  - )D  u]  - u = 0 (5.111) 

m^  + kY 

Letting  u = Dcj) , integrating  (5.111)  and  adjusting  the  in- 
tegration constant,  we  obtain 


D^(t)  + ^D4>  - (!^^  + k^)<)>  = 0 


<))  = A*I^(kr)  + B*K^(kr) 


(5.112) 


The  solutions  for  the  perturbation  velocity  and  pressure 
are 


^ = AkP(kr)  + Bkr  (kr) 
p = -ip^N^[AI^(kr)  + BK^(kr)] 


or  in  the  individual  region  we  have 


= A,  kr  (kr) 
N-|  1 m' 

I 

. P^  = -ip^N^LA^ 
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■< 


P2  = -iP2N^[A2l^(kr)  + B^K^Ckr)] 


(5.113) 


P3  ^ 


After  matching  both  the  kinematic  and  dynamic  boundary 
conditions  at  the  interfaces,  the  secular  relation  for 
stability  yields 


This  expression  can  be  reduced  to  that  of  the  two- region 
flow  by  letting  either  p.|  = P2  and  W-|  = W2,  or  03  = P3 
and  W2  = W3. 

Considering  p-|  = p^*  = W^,  i.e.,  an  annular 

jet  injected  into  a solid  jet  of  different  density,  we  ex- 
pand (5.114)  according  to  the  asymptotic  forms  of  the  modi- 
fied Bessel  functions  and  obtain  an  algebraic  relation 


(5.114) 
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[a(W2  - c)^  - (W^  - c)^]  - e^'''^[a(W2  - c)^  + (W^  - c)^]  = 0 


(5.115) 


Equation  (5.115)  can  be  solved  for  c 


c a 


( 1 + e*^^  ) T - (1  - e 


KE 


) 1 i ( 1 - T ) - 1) 


a(l  + e^")  - (1  - e^") 


(5.116a) 


(1 


KE 


)t  - (1  + e^^)  + id  - t)  Ja{e^^^  - 1) 


(1  - e^^)  - (1  + e^^) 


a 


(5.116b) 


Wo  W.  . , 

where  t = tA  = rA  , e = |-  , and  e is  not  necessarily  small. 

W3  R 

Equations  (5. 116a, b)  represent  unstable  solutions  existing 
at  the  interfaces  r = R and  r = R + 6 respectively.  For 
sufficiently  short  wave  lengths,  instabilities  will  take 
place  at  both  the  interfaces  for  t 1 , and  they  are  not 
expected  to  couple  with  each  other.  It  is  interesting  to 
note  again  all  the  unstable  solutions  lie  exactly  on  a semi 
circle  with  the  reason  given  earlier  in  the  discussion  of 
non-rotating  jets  of  the  first  example. 

Figure  5.26  shows  the  maximum  growth  rates  c-  = 

^|W2  - W-|  I in  the  a - E plane.  The  curves  in  the  lower 
part  of  the  plane  represent  the  solutions  for  maximum 
growth  rate  at  interface  r = R,  while  the  curves  in  the 
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Figure  5.26  Curves  of  maximum  growth  rate  (c^  = - W^|/2) 
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upper  part  of  the  plane  represent  those  at  the  interface 
r = R + 6.  They  both  approach  the  asymptote  a = 1 which 
is  the  most  unstable  situation.  Both  the  small  or  large 
density  ratios  have  stabilizing  effects  and  the  flow  will 
be  stable  for  a = 0 or  a = °°,  i.e.,  an  annular  liquid  jet 
directed  into  vacuum  or  an  annular  air  jet  blowing  into  a 
heavy  fluid  is  of  course  stable  in  the  absence  of  surface 
tension. 

5.4.2  Continuous  Axial  Jets 

The  second  example  of  jet  instability  we  wish  to 
discuss  is  a three-region  flow  with  continuous  velocity 
distributions  subject  to  axisymmetric  disturbances. 

Squire  (1933)  (see  also  Yih  1955)  has  proven  that 
plane  parallel  flows  are  most  unstable  to  two-dimensional 
disturbances.  It  is  suggested,  but  not  proven,  by  many 
analytical  and  experimental  results  that  the  two-dimen- 
sional azimuthally  periodic  perturbati ons  are  amplified  more 
strongly  than  the  corresponding  three-dimensional  ones. 

In  fact,  the  phase  velocities  of  the  perturbati ons  are  in 
general  expected  to  follow  the  direction  of  the  mean  flow. 
For  this  reason,  the  onset  of  instability  for  axisymmetric 
perturbations  would  be  expected  to  take  place  in  the  axi- 
symmetric jet  case,  and  should  therefore  be  discussed. 
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Consider  a velocity  and  density  profile  of  the 
following  form 

_ r 2 - ^ - 

W = A(^)  + B 

(5.117) 

rill'' 

Po  P2(r' 

here  P2,  and  a are  arbitrary  real  constants.  For 

axisymmetric  disturbances  ( m =0),  the  governing  equa- 
tion for  profile  (5.117)  taken  from  equation  (5.110)  is 

D(r'^D*u)  - k^r'^u  = 0 (5.118) 

a 

~2 

We  normalize  the  equation  by  setting  u = r ())  and  obtain 

(1  - -)^ 

d2^  + Id*  - [ ^ + k^]*  = 0 (5.H9) 

with  solutions  in  terms  of  the  modified  Bessel  functions 
of  the  first  and  second  kinds  of  order  y,  namely 

(})  = AI  (kr)  + BK  (kr) 

^ y y 


where 


y = n - fl 
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The  perturbation  velocity  is 

a 

u = r ^[Al^(kr)  + BK^(kr)]  (5.120) 

Following  Michalke  and  Schade  (1963),  we  construct 
a three-region  flow  with  velocity  and  density  distribution 
as  foil ows ; 


rW 


1 


W 


R2  ^ - r" 

- Rf 


+ W, 


0 ^ r £ R-| 
R^  ^ r < R2 


W 


3 


R2  1 r 


< (X) 
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0 < r < R-| 
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0 


= < 
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^3 
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r <00 
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r) 
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5.27 
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and  density 

distributions 
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After  taking  care  of  the  irregularities  of  the  modified 
Bessel  functions,  the  solutions  for  the  fluctuation  veloc- 
ities and  pressures  in  each  region  are 


< 


Pi 


P2 


A^kl' (kr) 

-iAi [p^N^ I^Ckr)] 


r ^[A2l^(kr)  + B^K^Ckr)] 


Pp  r ^ 

i4(f-)  r 2 


u 1 

a -75-  - 1 


(2  - a)  (W^  - W-|  ) 2 _ 


kr I ' ( kr ) 


(2  - a)(W2  - W^) 


I,TkTT^^^C"')  ^2^''n2  - a _ .2  - 


(5.121 ) 

^ - N[{1  - |)+ 

2 - a 
r 


krK'  ( kr ) 

y 

,U3  = B3kK;(kr) 

. B3[p3N3K^(  kr)  ] 


Matching  the  interfacial  condition  at  the  interfaces,  one 
finds  the  secular  relation  governing  stability  is 
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[F(k^)(W^  - c)  + ■ ^H]I^(k^)  [G(k^)(W^  - c)  + " ''h]K^(k^) 


2 - a, 


[F(k,){W,  - c)  + R2  ■ [G(k2)(W2  - c)  + ^H]K^(<2) 


'2' '"2 


where 


(5.122) 


F(k^) 


F(k2) 


G(<i  ) = 


G(k2) 


"1 

1 

0 

1— 1 

“1 

IC2 

Ko(i<2) 

a2 

Kjyqr 

^1 

“1 

Ko(i<2) 

w 1 / 

- (1 


- (1 


’ ■ "W 


I ( KoT 

y ' 2 


(1  - i-) 


'0^2 


- (1 


-) 

2' 


—lyTTy 


K2K^(i<2) 


H = 


(2  - o)(W2  - W^) 
R2  ■ ^ - R^  ■ ^ 


a- 


P O •] 
^1  *^2 


a. 


P3 


Equation  (5.122)  is  a quadratic  equati on  for  the  eigen  wave 
speed  c 


Ac 


2 


Be  + D = 0 
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or 

B + - 4AD 

c = 

2A 


(5.123) 


where 

A = F(k^)6(k2)I^(ki)K^(k2)  ‘ ^ ^1  ^ 

B = {F(k^)[W^G(k2)  + HR^  ■ + G(k2)[W2F(<i)  HRi  ” 

- {F(k2)[W^G(k^)  + HR^  ■ + G{k-|)[W2F(k2)  + HR2^  ^ 

D = [W.jF(ki)  + HR^  ^][W2G(<2)  + ^^2 

- [W^G(k^)  + Hr2  ■ ^][W2F(k2)  + HR^  " ^ 

Unstable  solutions  exist  whenever 

- 4AD  < 0 

This  condition  reduces  to  Michalke  and  Schade's  (1963) 
result  for  homogeneous  fluid  (a  = o)  and  W2  = 0.  They 
obtained  some  special  unstable  solutions  for  long  axial 
wave  lengths. 

If  we  expand  the  modified  Bessel  functions  accord- 
ing to  their  asymptotic  forms  and  assume  wave  lengths  are 
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sufficiently  small,  such  that  the  terms  involving 
g-2(ic2  - x-j)  negligible,  then  equation  (5.123)  can 

be  simplified  into 


F(k2)6(k^)c^  - {F(i<2)[W.|G(k-|)  + HR^  + G{k.|  ) [W2F(k2)  + HR2  ]>  c 


[W^G{k.|)  + HR^  ■ '^][W2F(k2)  + HR^  ' "^]  = 0 (5.124) 


In  this  case,  the  terms  inside  the  radical  are  given  as 

G(k-,)  F(<2)  2 

- 4AD  = (W,  - W„)  {F(k„)G(k,  ) + (2  - a)[  R,^  “ ° ~r7^^  - ° 

1 - (^)  ^ 

K2 

(5.125) 


Therefore  the  flow  will  always  be  stable  for  sufficiently 
short  waves.  This  result  is  expected  since  the  axial 
velocity  under  consideration  is  continuous  across  both  in- 
terfaces, and  therefore,  the  flow  is  stable  for  large  axial 
wave  numbers.  The  density  variation  in  this  case  does  not 
have  significant  influence  on  the  stability  characteristics 
because  of  neglect  of  gravitational  effects. 


CHAPTER  6 


CONCLUSION 

As  shown  in  the  preceding  chapters  both  via  some 
general  stability  criteria  and  by  exact  solutions  for  sever- 
al different  velocity  and  density  profiles,  the  instability 
mechanisms  present  in  heterogeneous  swirling  flows  may  be 
of  either  centrifugal  or  shear  origin.  The  classical 
Rayl ei gh-Synge  criterion  is  found  to  be  neither  a necessary 
nor  a sufficient  condition  for  the  stability  of  such  flows 
when  axial  shearing  rates  and/or  non-axi symmetri c disturbances 
are  considered.  A general  criterion  based  on  a non-dimen- 
sional local  number,  reminiscent  of  the  Richardson  number 
encountered  in  stratified  parallel  flows,  has  been  derived. 

It  is  a sufficient  condition  for  stability,  subject  to 
infinitesimal  arbitrary  perturbations,  and  is  found  to  never 
be  violated  for  the  specific  examples  considered  in  Chapter 
5.  Such  a sufficiency  condition  depends  completely  on  the 
existing  stabilizing  effects,  which  appear  mathematically 
as  a second  order  singularity  of  the  differential  equation 
governing  stability.  The  vanishing  of  such  a singularity 
(i.e.,  the  vanishing  of  the  stabilizing  factors  existing  in 
the  flow)  results,  automatically,  in  the  violation  of  the 
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sufficiency  condition  and  no  stabilities  can  then  be 
predicted  (the  same  argument  can  also  be  applied  to  the 
Tayl or-Gol dstei n equation  for  two-dimensional  stratified 
shear  flows).  For  this  reason,  no  stabilities  can  be 
guaranteed  for  flows  governed  by  stability  equations  with- 
out a second  order  singularity,  such  as  non-rotating  jets 
or  homogeneous  two-dimensional  parallel  flows  (i.e.,  Ray- 
leigh's equation),  before  solving  the  governing  stability 

equations. 

For  possible  instabilities  in  heterogeneous  swirling 
flows,  the  semi-circle  theorem  has  been  verified  to  be 
valid  for  four  types  of  shearing  flows  possessing  centri- 
fugally  stable  profiles.  Unstable  solutions  for  non-rotating 
jets,  rotating  flows  with  axial  velocity  gradients  subject 
to  axisymmetric  disturbances,  rotating  flows  for  azimuthally 
periodic  disturbances  in  addition  to  DCp^^^)  - 0,  and 
swirling  flows  for  the  narrow  gap  approximation  must  lie 
inside  a certain  semi-circle  in  the  complex  phase  velocity 
plane.  The  range  between  maximum  and  minimum  velocity  of 
the  steady-state  flow  serves  as  the  semi-circle  diameter. 
Another  type  of  semi-circle  bound,  proved  only  for  rigid 
body  rotation  subject  to  two-dimensional  azimuthally  peri- 
odic perturbations,  states  that  all  unstable  solutions  must 
lie  exactly  on  a semi-circle  in  the  complex  phase  velocity 
plane,  with  the  constant  angular  velocity  as  diameter. 

Several  exact  analytical  solutions  have  been  studied  and 
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the  results  have  been  brought  up  and  compared  with  predic- 
tions by  the  semi-circle  theorem.  It  is  interesting  to  note 
that  the  unstable  solution  for  a potential  vortex  surround- 
ing a non-rotating  jet  with  smoothly  varying  density  dis- 
tribution and  subject  to  arbitrary  perturbations  also  lies 
exactly  on  a certain  semi -circle  in  the  complex  c plane. 

The  velocity  and  density  profilesfor  this  case  however  do 
not  fall  into  any  of  the  categories  of  semi-circle  bound 
mentioned  above.  We  believe  that  the  flow  belongs  basically 
to  the  non-rotating  jet  type,  but  with  a superposed  potential 
vortex  which  has  stabilizing  effects  on  the  instability. 

The  reason  for  an  exact  semi-circle  in  this  case  is  mainly 
due  to  the  constant  axial  velocity  distribution  inside  each 
individual  region  as  mentioned  before. 

It  is  suggested,  but  not  proven,  by  the  analysis 
and  results  that  rotating  flows  with  discontinuous  velocity 
profiles  are  less  stable  with  respect  to  two-dimensional 
azimuthally  periodic  disturbances  than  to  the  corresponding 
three-dimensional  ones.  Shear  effects  in  this  case  play  a 
dominant  role  in  the  onset  of  instability.  The  corresponding 
most  critical  perturbations  for  non-rotating  jet  flows  are 
axi symmetric  ones. 

For  broken-line  profiles  with  velocity  or  density 
discontinuities  at  the  interfaces,  the  centrifugal  effects 
must  always  be  included  in  the  dynamic  interfacial  conditions. 
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Failure  to  do  this  led  Michalke  and  Timme  (1967)  to  perform 
an  erroneous  calculation  at  the  cylindrical  vortex  sheet. 

The  following  briefly  summarize  the  results  on  the 
stability  behavior  of  the  heterogeneous  swirling  flows  ob- 
tained in  this  study: 

1.  There  are  two  types  of  instability  mechanisms  for  hetero- 
geneous swirling  flows:  centrifugal  and  shear.  The  Ray- 
leigh-Synge  criterion  is  only  a condition  for  centrifu- 
gal stability. 

2.  The  sufficient  condition  for  stability  is  ^ [refer 
to  equation  (3.19)  for  definition]  everywhere  inside  the 
flow  domain. 

3.  The  semi-circle  theorem  holds  for  four  types  of  centri- 
fugally  stable  profiles.  Unstable  solutions  for  rigid 
body  rotation  against  azimuthally  periodic  perturbations 
must  lie  exactly  on  a semi-circle  in  the  complex  eigen 
velocity  pi ane . 

4.  Densities  increasing  radially  outwards  have  stabilizing 
effects.  This  has  been  noted  in  a related  study  by 

Ku rzweg  ( 1 969) . 

5.  Axial  velocity  gradients,  of  eigher  sign,  destabilize 
the  flow.  This  is  a manifestation  of  the  familiar 
Kel vi n-Hel mhol tz  instability  mechanism. 

6.  Azimuthally  periodic  disturbances  are  most  unstable  for 
rotating  shearing  flows. 
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7.  The  centrifugal  effects  due  to  the  discontinuity  of 
density  or  azimuthal  velocity  must  be  included  in  the 
dynamic  interfacial  condition. 

8.  Surface  tension  always  stabilizes  the  flow  except  for 
m = 0 when  k < 1 . 

9.  The  phase  velocities  of  perturbations  are  in  general 
expected  to  follow  the  direction  of  the  mean  flow,  with 
the  larger  wave  numbers  corresponding  to  stronger  ampli- 
fied growths  for  vortex  sheets  existing  in  the  flow. 

10.  Unstable  solutions  for  rigid  body  rotation  with  constant 

axial  flow  and  density  distribution  p^ir)  subject 

to  arbitrary  perturbations  are  solved  exactly  in  terms 
of  the  modified  Bessel  functions  with  complex  auguments 
and  orders.  For  azimuthally  periodic  perturbations  in- 
stability occurs  when  a < 0,  while  for  axisymmetric  per- 
turbations instability  will  not  take  place  until  a<-4. 

2 cj 

11.  Non-rotating  jet  flows  with  axial  velocity  W(r)  = Ar 
+ B and  density  PQ(r)~r'^  are  also  solved  in  terms  of 
the  modified  Bessel  functions  with  real  arguments  and 
orders.  This  profile  is  drawn  between  an  inner  and  an 
outer  concentric  jet  of  constant  density.  The  flow  is 
found  to  be  always  stable  for  sufficiently  short  axial 
wave  lengths. 
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